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1. STATEMENT OF WORK 

The statement of work for this project includes data analysis and 
supporting research in connection with the following broad objectives: 

(1) Provide a precise and accurate geometric description of 
the earth's surface. 

(2) Provide a precise and accurate mathematical description 
of the earth’s gravitational field. 

(3) Determine time variations of the geometry of the ocean 
surface, the solid earth, the gravity field, and other geo- 
physical parameters. 
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2. ACCOMPLISHMENTS DURING THE REPORT PERIOD 

2.1 Adjustment of the BC-4 Worldwide Geometric 
Satellite Triangulation Net 

2.11 Theoretical Developments 

, As was mentioned in the last semi-annual report, work was begun on 
processing the NOAA TYPE II data. The only existing computer program at 
that time was designed to use non-correlated data. The TYPE II data, being the 
result of a polynomial fit to plate images, has an associated 14 x 14 variance- 
covariance matrix, and in order to use this data it was necessary to write a new 
program. 

The new linear form of the mathematical model is 

F i . = Xg - X G - R cos a cos 6 

F 2 = Y* - Y 0 - R sin a cos 5 (1) 

F 3 = Z 3 - Z G - R cos 6 


where the subscripts S and G refer to satellite and ground, respectively, and R is 
the range from the ground station to the satellite. The observations are a and ff . 

The linearized form Of the mathematical model is basically the same as 
described in The Ohio State University, Department of Geodetic Science ReportNo, 86, 


(pp. 21-27), which is 


AX + BV + W = 0, 



where the matrices A and B are the partial derivatives with respect to the parameters 
and the observations, respectively. Whenever a satellite event is defined as the 
observations to one satellite position, the A matrix for one ground station and 
one satellite position is of the form 


+1 0 0 i -1 0 0 
0 +10 i 0 -1 0 

0 0 +1 i 0 0-1 



I] 


(3) 


A 


However, in case of correlated observations an event is defined as all obser- 
vations to the seven (7) satellite positions, and the A matrix for one ground 
station and seven satellite positions takes the form 



(21 x 24) 

This is perhaps easier to understand if the linearized form of the mathematical 
model is split up as follows: 

AiX a + AgXs + BV + W = 0. (5) 

This is essentially what was done in the original adjustment program. But 
when the model in the original program is split up, the A matrices are either 
+1 or -I and they cancel out in the mathematical development, the only change 
being that of signs. For the correlated data, Aj is the left side of equation (3), 
and it will not cancel out. 

Another change that had to be made was in the formation of the matrix 

M' 1 = (BP*B')*. (6) 

The problem arises because BP^B* is a singular matrix and cannot be inverted; 
For the case of one ground station and one satellite position one can use the 
following 

M' 1 = (BP rl BV = (B') 1 P(B) 1 = (B'VpB 1 , (7) 


where 




As can be seen in equation (7), the matrix B must be inverted, which means that 
it must be square. Thus even though the range R in equation (1) is not 
an observed quantity, it must be considered as such in order to 
make B a square matrix. This is of course compensated by inserting zeroes 
in equation (9). 

The above development for M' 1 is described in the above mentioned Report. 
In .case of correlated images the situation is somewhat more complicated. 

The matrix B is now of dimensions 21 x 21 and of the form 



where each of the blocks is a 3 x 3 as defined in equation (8). 

The matrix P cannot be defined quite as simply as in equation (9). The 
original variance -covariance matrix is 14 x 14, and the P matrix is 21 x 21.: 
This is handled as follows: 
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(21 x 21) 

By using the matrices B from (8) and P from (12), equation (7) can be 
solved for M’ 1 (the notation M' 1 is a misnomer, but this expression was used in 
Report No. 86 and we have continued with the same notation). The complete 
description fo the mathematical will be given at a later date. 
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By using the techniques described abouve, the reduced normal equations 
are formed as described in Report No. 86. 

In addition to the generalized approach described above, a completely dif- 
ferent mathematical model has also been developed using the method of observation 
equations. The principal advantage of the method of observation equations is that 
here the original given correlation matrix is used without any modifications which 
is necessary in the generalized least squares approach. 

2.12 Data Acquisition 

As of the end of this reporting period the following BC-4 data has been 

received from the data center: 

. / . • • . 

(i) Type I Data - 31 Tapes. 

(ii) Type II Data - 15 Tapes. 

The tape-wise details for type II data are as listed below: 


Tape Serial 
No. 


A-10806 

A-10268 

A-11082 

A-03725 

A-03719 

A-03727 

A-03728 

A-10897 

A-03738 

A-95575 

A-11519 

A-12327 

A-12037 

A-12010 

A-14094 


No. of events 
on the tape 


Break up of events with 
simultaneously observing stations 


2 stations 3 stations 4 stations 


















2.2 Investigations Belated to the Problem of Improving Existing 
Triangulation Systems by Means of Satellite Super-Control Points 

2.21 Introduction 

Geodetic triangulation has been accepted as an accurate method of 

determining "precise" coordinates for the triangulation stations of relatively 

short chains. This well-accepted idea was also given in an article "How 

accurate is First-Order Triangulation?" FSimmons, 1950, pp. 53-561 with 

the following introductory words: 

The question is often asked, "How accurate is the position 
of a triangulation station, " or "To what accuracy are the 
distances between triangulation stations known?" These 
questions are difficult to answer, principally because 
first-order triangulation gives the optimum accuracy 
in the measurement of great distances and there is at 
present no super yardstick to which it can be compared. 

Two modern technological advancements, namely, satellites and 
electronic distance measuring (EDM) instruments* have questioned the 
first-order triangulation accuracy, especially if triangulation is extended 
to distances longer than 1000 km or more. In such extended triangu- 
lation systems systematic errors like lateral refraction, propagation of 
observational errors, residual polar motion effects on latitude, longitude 
and azimuth, etc. [Mueller, 1969, pp. 61, 86-87; Pellinen, 1970, pp. 34-35; 
Wolf, 1950, pp. 117], which cannot be eliminated, accumulate. Lately 
the question has been raised whether any significant increment to accuracy 
is "cascaded" from a 1:1 million 1000 km net through a 100 km net to 
local control over 10 km distances. 

The satellite triangulation and super-transcontinental traverse, being 
of the highest achievable accuracy of today,, i.e. , super-control net of 
"zeroth" order, constitute a modern geodetic super structure, within 
which the classical geodetic triangulation is supposed to provide a geodetic 
control densification. 
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According to the classical geodetic ooncept, a lower order system 
should be lied to a higher ordor system. Statistically, this means that 
the variance-covariance of the higher order system, as a lower limit 
for accuracy, be at least compatible with the internal precision of the 
lower order system. For all practical reasons, the accuracy of the 
higher order systems should be substantially better (by a factor of two 
to three) than the subordinated system, thus supplying a rigorous con- 
straint in the reduction of the lower order system [Schmid, 1960, p. 4]. 


The objective of this investigation is to answer the question: 
Whothpr any significant increment to accuracy could be transferred 
from a supor-contvol net to the basic geodetic net (first-order triangu- 
lation). This objective was accomplished by evaluating the positional 
accuracy improvement for station Wyola (96), whloh is near the middle 
of the investigated goodotlc triangulation not, by using various station 


constraints over its geodetic position. 


2, 22 Data and Accuracy Estimates 

For the purpose of the present investigation, the triangulation of 
the western-half of the United States has been considered, as this is 
more accurate than that of the eastern-half of the United States 
[Simmons, 1950 , p. 54]. The investigation is done on the chain from 
Moses Lake, Washington to Chandler, Minnesote (Figure 1), .as these 
two stations are also common on both the continental satellite net (CSN) 
and the super-transcontinental traverse (STT). The data used were 
supplied by the Triangulation Branch of Geodesy Division, and the 
Geodetic Research ajid Development Laboratory, both of the National 
Oceanic and Atmospheric Administration, Washington. 






The details of Moses Lake - Chandler triangulation chain are as 
follows: 

Number of stations 101 

Num her of bases T Taped 27 

LGeodimetcr 2 

Laplace stations 13 

Observed directions 919. 

Distance between two stations [Minimum 273m 

t L Maximum 190 km 

Total length of the chain 1858km. 

It is assumed that the necessary reductions have been applied to 
the observed data, and the weight function P is "a priori" known to be 
a sufficient good accuracy. 

Super-transcontinental traverse (STT) runs across the western- 
half and the eastern-half of the U.S.A. (Figure 2). Its specifications, 
configuration, reduction of data and instrumentation are dealt with by 
Meade [1967; 1969a; 1969b}. 

Continental satellite net (CSN) is, in general, planned in such a 
way so that the stations are around 1200 km apart and that these stations 
are evenly distributed over the entire area. CSN-stations are either 
the stations of first-order triangulation net or these are connected to 
them. Its specification and configuration are dealt with in [Deker, 1967; 
Mueller, 1964; Pellinen, 1970; Schmid, 1970]. The continental satellite 
net of the North American Continent comprises of twenty stations which 
can be anchored in the three world net stations; Thule, Greenland, 

Moses Lake, Washington, and Beltsville, Maryland. Furthermore, 
planned is a tie to a fourth world net station - Shemya (Figure 3) 
[Schmid, 1970]. 

The following representative standard errors for observed data 
of Moses Lake-Chandler triangulation chain has been suggested [Meade, 
1970]: 
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Directions 

Azimuth 


0.4 

0."8 


Base 


Taped 

L Geodimeter 


1 pai l in f>00, 000 

l ppm for distance > 15 km 

1 . 5 cm for distances up to 15 km 


The mean of all section closures, which is the accuracy measure 
for the investigated geodetic triangulation net, is given as 1 part in 
317,000 [Adams, 1930], The standard position errors of the end 
stations of super-transcontinental traverse, which represent its accuracy 
measure, using actual data sets as given by different investigators 
differ too much from each other. The proportional error, which is the 
standard position error divided by the distance of the station from 
traverse-origin, is used for this investigation. The proportional 
errors of super-transcontinental traverse are given as follows: 

1:740,000 over 318 kilometer long traverse, and 1:1,100,000 over 1270 
kilometer long traverse [Foreman, 1970]; 1:670,000 over 270 kilometer 
long traverse [Gergen, 1970] and 1:3,000,000 over 1858 kilometer long 
traverse [ESSA, 1969]. The preliminary accuracy (i.e. proportional 
error) of continental satellite net, as obtained from the supplied data, 
corresponds to 1:385,000 for Chandler station. Because of this wide 
range in preliminary accuracy measures of these two super-control nets, 
investigations using the following accuracies (station constraints), are 
made: 1:300,000; 1:400,000; 1:500,000; 1:600,000; 1:700,000; 1:1 M; 

1:1.5 M; 1:3 M. The use of these accuracy measures, which are 
within the limits of preliminary accuracies of the two super- control 
nets, will determine a limit on the accuracy requirement of the super- 
control net, which would be necessary to improve the geodetic triangu- 
lation net. 
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2.23 C omputations and Results 

During the earlier period of this investigation considerable thought 
was given to the selection and use of such formulas and methods which 
would not only provide high accuracies, but also minimize or eliminate 
loss of accuracy in computations. This resulted in using Helmert- 
Itainsford-Sodano's Iterative Solution for Inverse Problem, which is 
equally applicable for short and long lines, and Conjugate Gradient 
Method (Cg- Method) for the adjustment of the triangulation nets, 
where the original observation equation coefficient Matrix (A -Matrix) 
is used, thus avoiding direct formation of normal equations where 
certain properties of the original A-Matrix are lost. To minimize 
the round-off errors, computations are done in double-precision with 
double precision storage [Muller- Merbach, 1970]. 

From the two basic adjustment methods, i.e. , Method of Obser- 
vation Equations and Method of Condition Equations, the former has 
been preferred for the present investigation due to reasons of simplicity 
and clarity. The reasoning of this preference has been dealth with in 
[Grossmann, 1961, p. 174; Helmert, I. Teil, 1880, p. 556; Wolf, 1968, 
p. 323] . Due to the large size of the triangulation net under investi- 
gation and the availability of digital computers, iterative methods were 
considered because (1) they are easier to program, (2) they require 
less storage space as the coefficient matrix of a triangulation net is 
very sparse, (3) they use directly the original set of equations through- 
out the process and hence rounding-off errors do not accumulate from 
one iterative cycle to another. 

While searching for a suitable adjustment method, this investigator 
came across the Conjugate Gradient Method (Cg-Method) [Schwarz, 1968 
and 1970; Wolf, 1968], which has the following distinct advantages over 
other iterative methods, such as Gauss-Seidl-, Jacobi-, Relaxtion- 
and other Gradient methods: 


- 16 - 



1. Original A-Matrix is used, thus avoiding the formation of normal 
equations, where certain useful characteristics of A-Matrix, such as 
very small coefficients may be lost. 

2 . Original A-Matrix, which has very few non-zero elements, is 
easily stored in comparatively much less computer space using an 
Index-Matrix. 

.V No "mesh-point numbering technique" TAshkenazi, 19671 to keep 
the band-width of the system a minimum is necessary. Thus stations 
can be added or taken out from the existing triangulation system with- 
out caring for their numbering. 

4. ft is a finite iterative process. Theoretically , the solution vector 
is obtained in a maximum of n-steps, n beingthe number of 
unknowns. Therefore, eigenvalues need not be calculated for 
determining the convergence. However, experimentation shows 
that the solution vector is not obtained in n-steps, as the 
orthogonality between the residue-vectors is not maintained 
exactly. Consequently, the residue-vector r tnJ after n-iterations 
is not zero. This deviation from zero depends upon the condi- 
tion of the system; the poorer the condition, the larger will be 
the deviation. 

6. Each approximation j$) to the solution vector is closer to the true 
solution x than the preceeding one. 

7. The ability to start anew at any point in the computation using the last 
x 0) as initial value so as^to minimize the effects of round-off errors. 

Following mathematical model; using method of observation equations, 
is used: 

Let Lj be the m independent observed quantities, v t the residuals to the 
observed quantities (obtained from the adjustment) and x,y, z, . . . the n unknown 
parameters to be determined. Each observation gives an observation equation, 
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whose general from is 

Lj + v, = f t (x, y, z, ... ), (1) 

where i = 1,2,3. . . , m and f represents a linear or non-linear function. The 
method of least squai’cs however demands that (1) f should be linear, i. e., 
a linear relationship between the observations and the unknowns and (2) the 
number of observations (m) should be greater than those of the unknowns (n) 
i. c.,m>n. In case of a non-linear function f this is linearized by using 
Taylor scries about such good approximate values of the unknowns x 0 , y 0 , z 0 , . . . 
such that the second and higher order terms can be neglected. In this 
case, equation (1) can be written as 

Vt = a 1 dx + btdy + c t dz + . . . + (2) 




It will be seen later that we have preferred to use weighted con- 
straints to the station Chandler. These "a priori" weighted constraints 
on the station position generate observation equations of the form 

v x = Fx (5) 

where F is a rectangular matrix, whose elements are either zeros or 
one. Thus the complete observation equation system can be written as 
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where 


V = Bx + L (6) 

v -u,j ! b •••[’!]• l = to j * 

Due to angular and linear (distance) observations, the observed data 
in a triangulation net are of a heterogeneous or dissimilar nature. 1 
This heterogeneous data have not only more than one dimension but also 
different "a priori" standard errors. To make this data homogeneous, 
i.e. , dimensionless and of unit weight, it is divided by the corresponding 
"a priori" standard error a. For reasons of simplicity, the mathematical 
model used is assumed to be uncorrelated. The resulting homogenized 
observation equation system can be written as 

V = Bx + L (8) 


where 


and 



v = v/cTt ; A = A/or t ; L .=? 1/Oj. 
v* = v x /a x ; F = F/cr x 


(10) 


0 ’ 1 = standard error of L t ; cr x = standard error of x y 


^’he term "heterogeneous or dissimilar" observations is used when the 
methods of their measurement are diverse; thus not only angles and dis- 
tances, but also distances and heights are heterogeneous observations 
[Wolf, 1968, p. 56]; [Schmid and Schmid, 1965a, p. 10] uses the term "hybrid 
systems" for "heterogeneous systems". 
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Equation (8) is used directly for adjustment by conjugate gradients method. 
A complete algorithm for obtaining solution vector and N* 1 by Cg-Method 
is given later, which gives v T Pv and Q*, or Q yy for a particular column. 
Using these quantities the "a posteriori" variance of unit weight (mf ), 
standard errors (m x , m y ) of unknowns, standard positional error (m p ) and 
the elements 0, A, B of the error ellipse are computed [Wolf, 1968, 
pp. 286-292]. 

\ The geodetic triangulation net is adjusted as an independent or free 
net, as it is not connected with other nets. For its unambiguous deter- 
mination, besides the observed data which include directions, bases (to 
provide the scale) and astronomical observations, i.e., longitude and 
azimuth (to provide orientation of the tir angulation net upon a mathe- 
matical surface, i.e. , ellipsoid), one fixed station is required to serve 
as origin [Gotthardt, 1968, p. 167; Grossmann, 1961, p. 175]. Moses 
Lake station is kept as origin with its coordinates obtained from satellite 
triangulation results; these coordinates have been assumed to be the best 
known coordinates. As Moses Lake station is fixed, its corresponding 
x-vector is zero, i.e., corrections dtp and dX are zero. For compu- 
tational ease their corresponding elements of the A-matrix are sub- 
stituted with zero. 

Combining the free triangulation net with super-control net of zero 
order, i.e. , continental satellite net and/or super transcontinental 
traverse means constraining the scale and/or orientation of the triangu- 
lation net. The effect of this combination is comparable with "tennis 
racket and string effect," where the rigid outer racket frame (super- 
control) constrains the loose strings (triangulation net). If the strings 
arc already constrained, there would be no "visible" effect of the 
additional constrain from the rigid outer frame. This is also the purpose 
of this investigation, i.e., to evaluate whether the existing geodetic 
triangulation is sufficiently "constrained" or needs to be constrained by 
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additional super-control net. For the present investigation triangulation 
station Chandler, which is common on the three networks, provides 
constraint. 

Geodetic triangulation net can be combined with the super-control 
net in either of the two ways : 

(1) By using the actual data, i.e. , by using the actual 
coox'dinates with their standard errors of Chandler as obtained 
from CSN and STT with the geodetic triangulation; or, 

(2) By adding a weight constraint to Chandler with its 
coordinates from the geodetic triangulation. 

For this investigation, the first way could not be used, as the super- 
control net coordinates of Chandler station are not compatible with those 
obtained from geodetic triangulation. As such, the second way has been 
preferred by using the actual preliminax'y accuracy estimates for Chandler, 
which are 1 part in 385,000 and 1 part in 3 million, as obtained from 
CSN ans STT, respectively. Further investigations are made by using- 
hypothetical standard positional error accuracy estimates of Chandler 
station, which are 1:400,000; 1:500,000; 1:600,000; .1:700,000; 1:1M; 

1:1,5 M. These accuracy estimates are within the actual preliminary 
accuracy estimates of super-control nets. Thus, using these various 
accuracies of super-control net, a feeling for the accuracy limit of super- 
control net, which would be necessary to improve the investigated geodetic 
triangulation, can be obtained. 

The Method of Conjugate Gradients (Cg-Method) is a nonstationary 
relaxation method, 

Nx + u = O (11) 

in n-iterative steps, where N is symmetric and positive definite. 

Then the system (11) - known in geodesy as the Normal Equations - 
has a unique solution. However, it is not necessary to have normal 
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equations, as Cg-Method can be easily modified for directly using the 
observation equations without explicit formation of normal equations. 

A complete mathematical derivation of Cg-Method with its program is 
given by Saxena [1972a; 1972b], 

A complete algorithm of Cg-Method for obtaining the solution 
vector (x) and for obtaining N" 1 using directly the homogenized obser- 
vation equations can be summarized in the following systematic way: 

A. For obtaining the solution vector (x) 

Given : Homogenized Observation Equation: Ax + 1 = v 


Select: Initial Trial Vector 

x l0J = 0 

Compute: 


(!) 

v (°) =» Ax( c ) +1 



Relaxation steps j - 1,2, 

n 

(2) 

r(i-») = A t v^) 


sov 

IAr^“ l ) T IAh( J “ l >) 

(for j * 2) 

(3) 

- (Ah^~^ T (Ah^ -l O 


, (i) 

1, ' ) •■L-r0- 1 ) . e 1 _ l h0- 1 ) 

(for j ~ 1) 
(for j * 2) 

P) 

J ■' (Ah( J )) r (Ah( j b 


m 

m 

xO) XjhO) 

y( J ) = Ax^) + 1 


(7a) 

v< cL>. :z v(J "' ) AhO) 


Tests: 



(8) 

Orthogonality Test: 



r(') T h( J ) = o 



r( J-l ) T r( J ). = 

0 
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( 9 ) 


r O~l)T 

003,51 ” ’ ‘ 

r ^ 1 ) 1 h^ 3 ) 

003 9a ” 1P55FP1 = 0 

v^) = 

cfrsck 


Termination of Iterations. 

Based upon the. theory of Cg-Method and the geodetic requirements, 
iterations should be terminated as soon as anj' of the following conditions are 
fulfilled: 

(a) if the improvement in the solution vector between two consecutive 

iterations is negligibily small, i. e., [x( J ) - - 1, 0-l(T 4 seconds 

(i. e. 1. 0*Kf 4 second in <p or X = 3.0mm), 

(b) if rtOMO « 0, 

(c) if (AhO)) T (Ah<’)) ■ *= 0; 

.(d) if the given number of iterations is reached; 

(o) if the round-off error (RFE) during iterations exceeds a certain 
accuracy limit, which is given by the vector difference 

l r t rue” r «J* where r t r , a A t Ax, +■ A T l = AM’) and r CoB? = A T v^ 6k - 

RFE> |A T (v(')- v (-0 k )| 3 

'Hie iterations should be terminated if r (J) 53 . RFE. 

B. For Obtaining N" 1 - Inverse of Normal Equations 

Given: Homogenized observation equation coefficient matrix A. 

Select: Initial trial vector = O; where q k is the k-th 

column vector of Q ( = N’ 1 ) 
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Compute: 

^ r ^ ~ “ e k J where e k is the k-th column vector of 
the unit matrix E . 

( 2 ) €j-i, h (J) , Xj are to be computed according to 
equations given in (A) above. 

(3) * X 3 h^ 

(4) r u } •=' r u_1> +XjA T (Ah u> ) 

Test and Termination of Iterations: Same as in (A) above. 

The algorithm of (A) is programmed as a SUBROUTINE SOLN and 
and (B) as a SUBROUTINE QSOLN. Both subroutines can be vised for 
any feasible size of data, which can be accommodated on the available 
computer, after changing KM, which is the PATSUM Basic Block Size 
for RTR. 

The main program used together with these subroutines has 
dimension statements and a data card for Number of Unknowns (NU), 
Number of Equations (NE) and Number of Columns of Index Matrix (NI), 
which can be changed if there is need for it. 

The program is universal in the sense that it can be used for 
varying data without much change and that "mesh-point numbering technique' 
is not required. Therefore, stations can be added or taken out from the 
triangulation system without worrying about the band-width and size of 
blocks. These programs have been tested on systems from as small 
as 2 unknowns, 3 equations up to as large as 804 unknowns, 1397 
equations. 

Although the Cg-Method theoretically gives the solution vector at 
n- iterative steps (n = number of unknowns), investigations show that the 
solution vector is not achieved in n-iterations due to round-off errors, 
ill-conditioning of the system, disturbances of the orthogonality and of 
the eonjugancy relations [Beckman, 1960, pp. 69; Hestenes and Stiefel, 

1952, pp. 411], The present investigation, using the actual data set, 
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shows that the number of iterations required to obtain the solution 
vector by Cg-Method using directly the A-matrix without explicitly 
forming the N-matrix depends upon two factors: (1) condition of the 

system, and (2) accuracy of the solution vector required. 

Using the geodetic triangulation data (573 unknowns, 963 equations), 
the program went up to 5778 iterations without giving any 7 decimal 
accurate solution vector, while 4 decimal accurate solution vector was 
obtained after 1161 iterations, i.e. , 2.1 times number of unknowns 
(Table 1). 

Each column vector q k of N" 1 is generally computed in less than 
1.2 n-iterations (Table 1). 

Table 1. 


Experiment 

Number* 

Number of 

Solution Vector 

Covariance Vector 
for Column 8 

Unknowns 

Equations 

Iterations 

Time** 
m sec 

Iterations 

Time** 
m sec 

1 

573 

963 

1161 

9 

37.13 

640 

3 

45.96 

2 

573 

965 

1177 

9 

23.27 

657 

3 

31.91 

3 

573 

965 

li75 

5 

45. 97 4 

659 

2 

12.59+ 

4 

573 

965 

1176 

9 

22.32 

682 

3 

45.64 

5 

573 

965 

1164 

5 

53.44+ 

.674 

2 

1.77+ 

6 

573 

965 

1162 

5 

41.16+ 

675 

2 

0.00 + 

7 

573 

965 

1166 

9 


631 

3 

20.03 

8 

573 

965 

1159 

9 

24.29 

648 

3 

ID. 29 

9 

J 

573 

965 

1169 

9 

29.41 

608 

3 

11.51 


+ Refer to Table 

♦♦Time is the Execution time on H-Compiler, Option = 2 (IBM 360/75) except 
those marked with a plus (+) sign, which is the Execution time on H-Compiler, 
Option " 0 (IBM 370/165). 
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The results of the investigation are given in Table 2 and 3, where- 
in the improvement of the particular geodetic triangulation by super- 
control net is visible only when its accuracy is at least 1 part in 500,000. 

Table 2. 


• 

Experiment 

Accuracy 

■ 


WYOLA (95) 



Number 

1 in 

m 0 

Q XX 

Qyy 


m y 2 

Remarks 

1 


2.42 

6.0 

0.5 

35.2 

2.9 

Free Net 

2 

300,000 

2.41 

6.7 

0.5 

38.9 

2.9 

. 

. 

3 

400,000 

2.41 

5.9 

0.5 

34.3 

2.9 


4 

’ 500,000 

2.41 

4. 1 

0.5 

23.8 

2.9 


5 

600,000 

2.41 

. 4.1 

0.5 

23.8 

2.9. 


6 

700,000 

2.41 

4.1 

0.5 

23.8 

2.9 


7 

1,000,000 

2.41 

3.7 

. 

0.5 

• 

21.5 

2.9 


8 

1,500,000 

2.41 

3.2 

0.5 

18.6 

2.9 


0 

3,000,000 

2.41 

2.1 

. 0.5 

12.2 

2.9 



Qyy and m,y ni y 3 are given in 10“ - seconds"'. 
















Tabic 3 


-M 

£ 



WYOLA (95) 

— 1 

Experime 

Number 

Accuracy 




Positional Improvement 
Relative to Experiment 1 

1 in 


m y 


Meters 

% 

B 

Free Net 

1.83 

0.37 

1.9 

. ■ 


B 

300,000 

1.93 

0.37 

2.0 

-0.1 

- 5 

9 

400,000 

1.81 

0.37 

1.8 

0.1 

5 

■ 

500,000 

1.51 

0.37 

1.5 

0.4 

21 

5 

000,000 

1.51 

0.37 

1.5 

0.4 

21 

6 

700,000 

1.51 

0.37 

1.5 

0.4 

21 

■ 

. 

1,000,000 

1.43 

0.37 

1.5 

0.4 

21 

B 

1,500,000 

1.33 

0.37 

1.4 

0.5 

26 

m 

3,000,000 

■ 

1.08 

0.37 

• 

1.1 

0.8 

42 


Standard Errors of Unknowns (m x , rn y ) and Standard Positional 
Error (m p ) are given in meters. 
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500,000 





Worth mentioning is that the longitude terms, which are Q yy and 
m y 3 in Table 2 remain practically uneffected. This could be explained 
by the fact that station Wyola is very close to Laplace stations, which 
control the azimuth error accumulation and thus control the longitude 
terms. 

It is interesting to note from Table 3 and Figure 4 of the inves- 
tigated geodetic triangulation chain lies between 1:300,000 and 1:400,000, 
which is quite in agreement with its presumed accuracy of 1:317,000. 

2.24 Summary and Conclusions 

The super-control net, i.e., continental satellite net or super- 
transcontinental traverse; can provide a useful constraint to the inves- 
tigated geodetic triangulation net, and thus can improve it only when the 
accuracy of super-control net is at least 1 part in 500,000; in this case, 
this corresponds to ±3.7 m standard position error for the station 
Chandler. 

The preliminary accuracy of super-transcontinental traverse is 
already better than this limiting accuracy of 1 part in 500,000. The 
preliminary accuracy of continental satellite net is, however, lower than 
the limiting accuracy of 1:500,000; the preliminary standard position 
error for Chandler as obtained from continental satellite net corresponds 
to ±4. 8 m, i.e., 1:385,000. The future will show whether the limiting 
accuracy could be achieved by continental satellite net, especially because 
numerous spatial triangulations of CSN have produced accuracies within 
the range of 1 part of 400,000 and 1 part in 700,000 [Schmid, 1965, p.22], 

Schmid [1970, pp. 23—24] indicates that continental satellite net will 
fall short on an optimum solution with respect to both its coverage and 
its accuracy. The three-dimensional positions of CSN-stations will 
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probably be determined to no better than ±4 meters in all components, 
which does not seem to be good enough at least for this particular 
investigation. 

ft might be useful to have a "block constrain" instead of "chain 
constrain", that is, to use four well separated satellite stations, namely 
003, 102, 112 and 134 (Figure 1). 

Super-transcontinental traverse can provide a better constraint, if 
more than two of its stations are common to the stations of geodetic 
triangulation net. Also, a ’block constrain", as explained above, might 
be more useful instead of a "chain constrain". 

The development tendencies of instrumentation indicates that the 
future super-control nets will use VLBI (Very Long Baseline Inter- 
ferometry) and Laser ranging systems. 
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2. 3 Geodetic Satellite Observations in North America (Solution NA9) 

The coordinates of several tracking stations tied to the NAD datum were 
computed through available observations to the GEOS-I satellite. Up to date 
the NA6 adjustment [Mueller, Reilly and Schwarz, 19691 and NA8 adjustment 
TMueller, and Reilly, 1971] had been published. The latter solution was performed 
using height constraints deduced from the SA069 geoid [Gaposehkin and Lambeck, 
19701. 

Recently a new detailed geoidal map with claimed accuracies of ±2 m, 

(on land), based on gravimetric and satellite data, was presented ("Vincent, Strange 
and Marsh, 1971]. With the new geoid, and the orthometric heights given in 
[NASA, 1971] more reliable height constraints were calculated as follows: 

From the initial values of the shifts SAO-NA8 (computed using the published 
shifts SAO-NAD and NAD-NA8 in [Mueller and Reilly, 1971]) and by an iterative 
process self-explained in Figure 1, the initial NA8 rectangular coordinates were 
shifted to the SAG origin and the geodetic coordinates computed. The ellipsoidal 
heights then were constrained using the undulations from [Vincent, Strange and 
Marsh, 1971] . With the original <p and A and this new height a new set of rec- 
tangular coordinates was obtained. Following this procedure iteratively, 
several shifts of this kind to the "geocenter" were performed until the sum of 
the undulation differences was very small. Through this process "best" shift 
to the geocenter was obtained. This shift was also used to compute the pre- 
liminary coordinates to obtain the reduced normal equations for the MOTS and 
PC-1000 optical data in the solution MPS7 ([Mueller and Whiting, 1972] and 
Section 2. 4). 

At all stations, a weighted height constraint was imposed, after shifting 
(with the above obtained values) to the final "geocentric (GC)" coordinates. Also, 
as in the NA8 adjustment, a distance constraint was imposed between stations 
3861 and 7043. Due to a recent correction in their coordinates a difference of 
3 m from the previously used value was taken into consideration [Meade, 19721. 
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Finally, unlike in the NA8 solution, "inner adjustment constraints" were 
also imposed in order to define the origin of the system in its most favorable 
position from the error propagation point of view fBlaha, 1971]. 

The coordinates of the NA9 solution are presented in Table 1 with cor- 
responding standard deviations. The coordinates transformed to the NAD datum 
are in Table 2. 

Table 3 shows the constrained heights at each station and the final undula- 
tions compared with those published by Vincent et al. In column AN and in 
parenthesis, the differences published previously in [Vincent, Strange and Marsh, 
1971] are shown. It can be seen from Table 3 that only two stations show sub- 
stantial disagreements (3903 Herdon, Virginia and 3407 Trinidad). It seems 
clear that the orthometric height as given in [NASA, 1971] for the station 3903 
has a gross error. Appropriately, the NASA Directory of Tracking Stations 
points out in the description of the referred station: "coordinates unverified, 
survey details are lacking. " The discrepancy with respect to station 3407 may 
be due to the fact that it is situated in the Caribbean, where large geoidal 
gradients are present. 

Table 4 shows the transformation parameters between the different systems. 
Listed on the first page are the 3 parameter-transformation solutions (only shifts 
considered), and the general 7 parameter solution. In this latter transformation 
the rotations were first computed through direction cosines independent of 
translations and scale factor (see Section 2.5). These rotation parameters con- 
strained with their variances were used in the final solution shown pp. 2-5 of 
Table 4 with the resulting variance-covariance matrix and the correlation co- 
efficient matrix for each transformation. In the variance-covariance matrix the 
angular units are in radians. 
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Figure 1 continued 
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Re: Table 3 


New Height 
Constraints 
h 


Re: Table 1 


Program 

"OSUGAP" 

Solution of 
Normal Equatj 

(X, Y, Z) NAD 


Program 
"DATMTR" 
NAD - NA9 


(X,Y,Z) NA9 

(X,Y,Z) nAD 


NA9 NAD 
(X.Y.Z) 


NA9 -»'NAD = 


Re: Table 2 



X 

Y 

.ZJ 


( a » k) N AD 


(X, Y, Z) NA9 
(X,Y,Z) gC 


NA9-»"GC" = 


Figure 1 continued 


AX = -6.7 m 
AY = 0.2 m 

AZ =■ 0.4 m 

NA9 + (NAD - NA9) 

-(p- 

X 

. h. 

AX = -51.7m 
AY = 144.1m 
AZ = 210.5 m 

: NA9 + ("GC" - NA9) 
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Table 1 

Coordinates of the North American GEOS-I Tracking Stations from the NA-9 Geometric Adjustment 
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All coordinates and standard deviations in meters. 
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General Information: No. of degrees of freedom 5228 

No. of ground stations 34 Quadratic sum of the residuals (V Vv) 5039 

No. of spatial chord equations 1 Standard deviation of unit weight 1 





NAD Coordinates of the North American GEOS -I Tracking Stations 
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Note : The above coordinates were arrived at by applying the shifts AX --6.7m, AY = 0. 2m and 
AZ = 0. 4m to the NA-9 coordinates and then converting these values to ellipsoidal coordinates on 
the ellipsoid a = 6378206. 4, b = 6356583. 8. 
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Table 3 

Height Constraints and Undulations 
(all units in meters) 


Number 

Station 

Constraints 

N 

AN 

h 

a 

"GC"* 

[Vincent et al. ] 

1021 

Blossom Pt. , Md. 

9 

3 

- 27 

-26 

- 1 

(" 7) 

1022 

Fort Myers, Florida 

23 

3 

- 16 

-18 

2 

< 1) 

1030 

Goldstone, Calif. 

898 

3 

- 23 

-27 

4 

( .8). 

1032 

St. John's, Nswf. 

102 

5 

12 

13 

- 1 


1033 

Fairbanks, Alaska 

165 


16 




1034 

E. Grand Forks, Minn. 

256 

3 

- 13 

-18 

5 

( 11) 

1042 

Rosman, N. C. 

916 

3 

- 23 

-22 

- 1 

(" 3) 

3106 

Antigua, W.I. 

8 

3 

- 45 

-40 

- 5 

(- 2) 

3334 

Stoneville, Mississippi 

45 

3 

- 20 

-19 

- 1 


3400 

Colorado Springs, Col. 

2184 

5 

- 4 

-10 

6 


3401 

Bedford, Mass. 

89 

5 

- 27 

-21 

- 6 


3402 

Semmes, Alabama 

84 

3 

- 21 

-18 

- 3 

(-12) 

3404 

Swan Island 

79 

7 

^ 32 


• ' 


3405 

Grand Turk, B. I. 

0 

5 

- 51 

-47 

- 4 

(-27) 

3406 

Curacao, N. Antilles 

44 

5 

- 30 

-26 

- 4 


3407 

Trinidad, Tobago 

285 

5 

- 50 

-34 

-16 


3648 

Hunter AFB, Georgia 

19 

3 

- 19 

-24 

5 

(- 5) 

3657 

Aberdeen, Maryland 

7 

3 

- 27 

-26 

- 1 

(- 4) 

3861 

Homestead, Florida 

16 

3 

- 22 

-22 

0 


3902 

Cheyenne, Wyoming 

1882 

5 

- 8 

-10 

2 


3903 

Herndon, Virginia 

132 

3 

-100 

-26 

-74 


5001 

Herndon, Virginia 

132 

3 

- 27 

-26 

- 1 


5333 

Stoneville, Mississippi 

45 

3 

- 17 

-19 

2 


5649 

Hunter AFB, Georgia 

23 

5 

- 23 

-23 

0 


5861 

Homestead, Florida 

22 

3 

- 27 

-22 

- 5 

<- 5) 

7036 

Edinburg, Texas 

72 

3 

- 8 

-11 

3 


7037 

Columbia, Missouri 

270 

3 

- 17 

-24 

7 

( 10) 

7039 

Bermuda 

26 

3 

- 38 

-36 

- 2 

(-2) 

7040 

San Juan, P. R. 

57 

5 

- 40 

-41 

1 

( 5) 

7043 

Greenbelt, Maryland 

56 

3 

- 30 

-26 

- 4 

(-15) 

7045 

Denver, Colorado 

1787 

3 


-13 

6 

(16) 

7072 

Jupiter, Florida 

26 

3 


-24 

1 

( 1) 

7075 

Sudburg, Canada 

276 

3 


-31 

3 

( 20) 

7076 

Kingston, Jamaica 

473 

3 


-23 

3 

( 20) 


*The geocentric coordinates were obtained from the NA-9 by adding the following 
shifts: AX = -51.7m, AY - 144.1m, AZ = 210.5m-. 
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Table 4 

Transformation Parameters 


cd 


TO 


U 

* 9 

a ^ 


No. Stations 
AX(m) 

AY(m) 

AZ(m) 

AX(m) 
AY(m) 
AZ(m) 
0x(") 
9y( // ) 

9. (") 

efxKJ 6 ) 


c 

Cd 

X-i 

U 

a > 

-u 


cd 

x-i 

cd 

a 


NA9-NAD 

32 

6.7 ±1.3 
- 0.2 ± 1.1 
-0.4 ±1.2 

-1.9 ±3.0 
-20.9 ±5.0 
23.5 ±4.3 

- 0. 80 ±0. 09 
0.56 ±0. 07 

- 0.25 ±0.11 

- 4. 82 ±0. 89 


NA9- n GC" 

■■'34''. 

51.7 ±0.7 

-144.1 ±0.8 
-210.5 ±0.9 

39.1 ±1.7 
-144.4 ±4.0 
-202.5 ±2.8 

- 0.37 ±0.04 

- 0.22 ±0.03 

- 0.22 ±0.05 

- 0.48 ±0.72 


NA9-SAO 

11 

35.3 ±2.6 

-148.3 ±2.6 
-175.0 ±2.6 

24.9 ±9.4 
-200. 0 ±11. 7 
-173.5 ±11.1 

- 0. 66 ± 0.29 
0. 14 ± 0.24 
0. 94 ± 0.35 

- 6. 78 ±1.88 


”GC"-NAD 

32 

-44.8 ±1.3 

144.8 ±1.1 

209.8 ±1.2 

-37'. 8\ ±2.6 
124.3 ±5.0 

227.7 ±4.0 

- 0. 35 ±0. 07 
0. 84 ±0. 06 

- 0.10 ±0.09 

- 4.28 ±0. 89 


♦Rotation parameters constrained (see Section 2.5) 
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T able 4 continued ■ ■ , . ■ ■ ■ ' 

ROTATION PARAMETER S CONSTRAINED 


r ~“V'GC"-NAD 

* ‘ “ **"£ -V.,; ' . ' *. 


SOLUTION FOR 3 TRANSLATION, 1 SCALE AND 3 ROTATION PARAMETERS 


_ax.._ JXL. A z . c ..... 0 -y- Ox- 

meters METERS METERS xlO**) SECONDS SECONDS SECONDS 

LEit-gJit J! 2i-gA6-_.. _._z A 2_g -O.JO 

__ - • . 4 ‘ , 7 VARIANCE - COVARIANCE MATRIX 

Q^jS^jL+QJ Q^JLaiQ±Q.Q QjLM^P<L..-P^3JiMr-Qb . - 9 • T5itJ>- 06 

0.T03D+0Q . 0.2460 + 02 -0.115D+02 0.40 9 D-Q5 0.242D- 06 „ . 0 . 9780- 07 -0.59 7..D- 06 

0.33 9D+Q1 - 0«1 15.D+.C 2 0...JjSJJT+.0 j 2__-0-,ZS 9 D-J 05.__JT-MlfcC6 0..LL1 D- 0.6 -0.9 6 0 D- 06 

zJ^33ilDr-Cl6 7.D.-rX2---Q-^.JLPnM 0^12-Q-ll ... 0.9 BO. p- 1^ 

0.6750- 06 0. 2420-06 0. 38^06- -0.83 1 D-LA 0.,l..L7ib 1 2„ Q.204D-1 3 - 0 . 7 0 2 p- 1 3 

0.3 5 8 0— 06 0 . 9 78 D— 07 0 . 1 1 lD^-06 0 . 3 12 0-15 0.204D-13 ..0 ...7;2.4D- 13 - C . 2 60 1> I 3 

Q-.-9_8J)D - 1. Q. 7 0.2 Pr 1 3. .^^2600- 1_ 3_ . 0 . 1 8.3 0- 12 


COEFFICIENTS OF C Q RRE LATION 



- . Q^01 0 r p2’ Q;.1 : 60DT 01 -Q.«ftQ6+0C£ 0.923 C + 00 

- 0*1^30+00 : 0.92 30 +0 0 -0. 8060+00 0.1 00 0+01 

, 0. 76^5D-*-00 0; i42'p -»-00 0. 278D +Qb --0 .272 C-01 

__Q,.3:L7D+ 0Q 0.7320-01 0. 1 0 30+00; C. 1300-02 


0.765D+ 0Q 0_.51.7D + 00 -0 .412.0 + 00 

0.1420+ 00 ;. .0 . 7 320-0 1 . -0.2310+00 
0 .278 0 + 00 0. 1 C 3 D + 00 ~0. 5 59.0+00 

-0. .2720-0.1 ' 0.130D-0 2 0.2 57 0-01 

0. 10 00+01 0. 2210+00 -0.4790+00 

0.2210+00 C .100D + Q1 - 0.2 26.D+ GO 

-0.4 790+00 -0.2260 + 00 . 0.1000 + 01 
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SOLUTION FOR 

3 TR ANSI. AT I ON, 1 SCALE AND 3 ROTATION PARAMETERS 
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• 

■ 
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AX ,AY AZ 

> e 

9. 

9„ 0v 



METERS MET 

ERS METERS 

; t> • • 

• ( xlO” 6 ) 

SECONDS SECONDS SECONDS 


24.89 -200 

.00 -173.46 

-6.70 

-0.66 0.14 0.< 

?4 










VARIANCE - 

COVARIANCE 

MATRIX 



0.887D+02 

O.165D+02 

0. 2300+02 

-0.152C-05 

0.11 ID- 04 

0.7C1D-05 

-0.515 D-05 

0^06512+122. 

Q .13.8.0 +0.3. 

— 0 . 1.1 IP + 02. 

0.1780-04 

0.394 D- 0 5 

0.3 3 4 D-05 

’-0.110 D- 04 

.0.2800+0.2 

-0.111.0+02. 0«1.24 0+03 

-0 .1 29C-04 

0.381 D-05 

0.850D— 06 

-0.142D-04 

-0. 1520-05 

0.178.0-04 

-0.1290-04 

0.354C-1I 

0.116D-13 

C. 766D-14 

-0.230D-13 

0,11X0-04. 

0.3940-05 

0.38 ID- 05 

0.1 160-13 

0.192D-11 

0.405D-12 

-0 . o 1 1 0- 12 

0*7010-05 

0.1 34 D-05 

0. 8500-06 

0.7660-14 

0.4050-12 

0.1350-11 

-0. 3040-12 


-0.1100-04 

-0. 1420—04 

-0.230 C— 13 

-0. 8110-12 

-0.304D— 12 

0.2H9D-11 

— : : 


COEFFICIENTS OF CORRELATION 
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0.1 00D+01 

0 . 1 490+00 

0.266D+00 

-0 .858C-01 

0. 848D+00 

0.64CD+00 

-0. 322D+00 

0. L4.9JI+.0.y. 

0 .1.0,0 0 +.01. 

-0. 84 70- 01 

>0.8060+00 

0.242D+00 

0,9800-01 

-0.5 530+00 

0* 2660+00 

- 0.3 47 D- 01 

0. 1000+01 

^0 .61 5C+00 

0 • 246D+00 

0. 6550-01 

— 0.7 480+00 

. -CL fl .8.5..8L0*jP.i. 

O.8Q60+00 

-0.615D+0Q 

0.1000+01 

0.444D-02 

0.350D-02 

-0.721 C— 02 

„a..8A8,D+.0.0. 

0.242.0+ 00,., 

0. 2460+00 

0 .444 C— 02 

0. 1000+01 

0. 2510 + 00 

-0 . 344Q+00 

Q.64QD+QQ 

0.9800-01 

0.655D-01 

0 .3500-02 

0. 2510+00 

0 . 1 OCO + Cl 

-0. 1 540+00 

-0^32.2.0 + 00. 

-0 . 5 53D+QO 

-0.7480+00 

-0. 7210-02 

-0 .3440+ 00 

-0. 154D + 00 

0.1 00 D + 0 i 
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ROTATION PARAMETERS C ONSTRA I NED 


- NAD . ' ’VA-: 


SOLUTION- FOR 3 TRANSLATION, 1 SCALE AND 3 ROTATION PARAMETERS 


>■>./% . . ; ..n-fe-’v • 

AX AY .. AZ ; 6 Z .. 0 y , _ 0 X 
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VAR. I ANC 6 r COVARIANCE 

MATRIX 
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• " • ' ’• • . ! . -‘;V. ' • ” **•.*• t 

0.1430-06 0.1630^06?. 0.447D-15 
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-0. 3600- 13 

-0. 665D-06 

— 0. 874D— 06 -0. 141CF-05* 0.144C-13 

-0. 1030-12 

— 0.380 D- 13 

0 .2680-12 


’■ ?.■ ■. . . y r". - • ■ ■ ■ •• • • 

7 . • 
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2. 4 Free Adjustment of a Geometric Global Satellite Network 

(Solution MPS7) 

2.41 Introduction 

The basic purpose of this experiment was to compute reduced 
normal equations from the observational data of several different 
systems described below to combine them eventually with the normal 
equations of the Wild BC-4 observations taken in the DOD/DOC 
cooperative worldwide geodetic satellite program and provide station 
coordinates from a single least squares adjustment. The solution 
described in this paper is a partial one obtained without the use of 
the BC-4 data. The observational systems combined were the 
Baker-Nunn simultaneous camera observations from the SAO world- 
wide network; the C-Band range observations from the NASA net- 
work; the MOTS and PC-1000 optical observations in North America; 
miscellaneous camera observations in Europe which were included 
in the SAO 69 solution; _ and, lastly, a group of optical observations 
where Baker-Nunn cameras observed simultaneously with MOTS 
and/or PC-1000 cameras in the previously mentioned group. 


2.42 Description 

Smithsonian Data 

A set of optical observations were obtained from the Smithsonian 
Astrophysical Observatory. These included 14,356 simultaneous obser- 
vations from 28 stations in the SAO 69 Network. For each observation 
the track angle was provided along with the uncertainties along and 
across the track. The variances and covariances, in terms of right 
ascension and declination, were computed as described in [Girnius and 
Joughin, 1968] . 

MOTS and PC-1000 Data 

The set of optical observations used here were the same as those 
used in the NA6 adjustment described in [Mueller, et al. , 1969]. The 
observations had been previously screened and a set of reduced normal 
equations, referenced to the North American Datum, obtained. 

In the meantime [Vincent, et al. , 1971] published a geoidal map 
based on gravimetric and satellite data. By an iterative procedure a 
new solution was computed which constrained the new undula- 
tions, and thus a set of geocentric coordinates were obtained. With 
these coordinates as initial values, but with the original set of 
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observations, new reduced normal equations were computed to be used 
in the solution described in this paper. 

' C-Band Observations 

The C-Band solution is a least squares adjustment of the range 
observations from twenty- eight C-Band radar stations operated by NASA 
in a worldwide network, which resulted in distances between the stations 
and a set of coordinates of the stations on the SAO C-G ellipsoid along 
with their standard deviations [Brooks and Leitao, 1970]. Upon request, 
NASA/Wallops Island kindly sent us the correlation matrix for these 
solutions, which enabled us to reconstruct the full variance-covariance 
matrix. 

Some of the stations in this adjustment could be related through 
ground triangulation to nearby Baker-Nunn, MOTS or PC-1000 cameras, 
thus the interstation distances provided indirectly the scale of the 
solution. The C-Band data was treated as though they were length 
observations between the stations and developed a program that com- 
puted the normal equations that would correspond to these length- 
observations utilizing also the reconstructed variance-covariance 
matrix. 

The computed lengths are listed in Table 1. 


Table 1 


Stations 

Length (m) 

Merritt Island (4082) to Pretoria (4050) 

10,909,592 

Merritt Island (4082) to Kauai (4742) 

7,362,142 

Merritt Island (4082) to Bermuda (4740) 

1,593,106 

Merritt Island (4082) to Grand Turk (4081) 

1,230,691 

Merritt Island (4082) to Antigua (40G1) 

2,288,026 

Kauai H. I. (4742) to Vandenberg AFB (4280) 

3,977,684 

•• 
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Mixed Optical Observations 

We also received from the NASA, National Space Science Data 
Center a magnetic tape containing records of optical observations on 
GEOS-I from November of 1965 to August of 1966. These included 
2322 simultaneous observations between Baker-Nunn cameras, MOTS 
cameras and PC-1000 cameras located on and around the North 
American continent. 

It being intended to combine the normal equations developed from 
the above observations with a set of normal equations developed from 
a much larger set of SAO, MOTS and PC-1000 described above, observations 
the possibility of duplicating observations had to be considered. 

In the case of the majority of the MOTS and PC-1000 and all of 
the Baker-Nunn observations , the few duplicated observations were 
overwhelmed by the large number of other observations at these 
stations. However, a number of MOTS and PC-1000 stations in the 
Caribbean area contributed only a few observations to the North 
American normal equations. Any duplicated observations here would 
have had an inordinate effect. Therefore, all such observations were 
eliminated. 


2.43 Constraints 

Inner Adjustment Constraints (Free Adjustment) 

The large number of optical observations effectively determined 
the orientation of the total network while the C-Band observations 
provided a scale. Only the origin remained undetermined. To 
define the origin of the system in its most favorable position (from 
the error propagation point of view) we imposed "Inner Adjustment 
Constraints " compelling the trace of the variance-covariance matrix 
to be a minimum [Blaha, 1971]. 

_ Length Constraints 

The C-Band observations described earlier introduced scale into 
our adjustments. They also provided much needed extra connections 
from Africa across the Atlantic and to the Caribbean Islands, and 
the length Kauai to Vandenberg Air Force Base greatly strengthened 
the geometry in the western United States 

In addition to the C-Band scale wc also introduced a weighted 
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chord length constraint between Homestead, Florida and Greenbelt, 
Maryland derived from updated Cape Canaveral datum coordinates of these 
two stations determined from the high precision geodimeter traverse 
in the eastern United States. 

Height Constraints 

At all stations, a weighted height constraint was imposed. The 
heights above mean sea level were obtained from [NASA, 1971] and 
to these, the undulations referred to the SA0 69 ellipsoid were added. 
The undulations were determined from a number of sources. Between 
North America and Europe the geoid of [Vincent, et al. , 1971] was 
used. In this report, the undulations of some stations were also 
tabulated (computed). These tabulated values were constrained with 
weights corresponding to a standard deviation of 3m. Other station 
undulations were intex'polated from the geoid map itself and, allowing 
for interpolation errors, received assigned standard deviations of 5 m 
except in those areas near the Caribbean where, because of large 
geoidal gradients, a standai’d deviation of 8m was estimated. For 
stations in other parts of the world (not covered by the above geoid 
map) the undulations were obtained from the SA069 geoid map, and 
standard deviations from 8 m to 15 m were assigned depending upon 
the number of gravity measurements available in the surrounding 
area. All heights constrained (H) are shown in Table 2. 

These height constraints, which are in effect independent obser- 
vations, provided a valuable strengthening of an otherwise weak 
geometric network. A test adjustment was run (MPS9) in which all 
previously described constraints were held except the height constraints 
and in this adjustment the final standard deviations of the coordinates 
were more than doubled and at poorly determined stations more than 
tripled. 

• • - Relative Position Constraints 

These weighted constraints were used to tie together the C-Band 
radar stations with nearby camera stations through the connecting 
triangulation, and also helped to connect the Baker-Nunn stations with 
nearby MOTS and/or PC- 1000 stations. 

In every case, Cartesian coordinate differences were computed 
on the local datum and the weights determined from standard 
deviations computed from a formula given in [Simmons, 1950]. 

This estimate was used in all cases except between Merritt Island 
and Jupiter, Florida, where the uncertainty was estimated to be 
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Table 2 


Height Constraints and Undulations 
(all units in meters) 


Number 

’ 

Station 

Constraints 

N 

11 

a 

MPS7 

F Vincent et al. , 19711 

1021 

Blossom Pt. , Md. 

- .20 

3 

- 30 

-26 

1022 

Fort Myers, Fla. 

- *3 

3 

- 18 

- 18 

1030 

Goldstone, Cal. 

902 

3 

- 27 

- 27 

1032 

St. John’s, Nswf. 

82 

5 

12 

14 

1033 

Fairbanks, Alaska 

188 

15 

4, 


1034 

E. Grand Forks, Minn. 

238 

3 

- 15 

- 18 

1042 

Rosman, N. C. 

887 

3 

- 24 

- 22 

3106 

Antigua, W.I. 

- 37 

3 

- 41 

- 39 

3334 

Stoneville, Miss. 

20 

5 

- 20 

- 19 

3400 

Colorado Springs, Col. 

2173 

5 

- 8 

■ 

- 11 

3401 

Bedford, Mass. 

63 

5 

- 28 

- 20 

3402 

Semmes, Alabama 

55 

3 

- 24 

- 18 

3404 

Swan Island 

31 

15 

- 38 


3405 

Grand Turk, B. I. 

- 29 

3 

- 39 

- 31 

3406 

Curacao, N. Antilles 

- 19 

8 

- 29 

- 26 

3407 

Trinidad, T. & T. 

221 

8 

- 59 

- 34 

3648 

Hunter AFB, Georgia 

- i2 

3 

- 25 

- 24 

3657 

Aberdeen, Md. 

- 20 

3 

- 26 

- 26 

3861 

Homestead, Fla. 

- 22 

3 

- 24 

- 22 

3902 

Cheyenne, Wyo. 

1872 

5 

- 11 

- 10 

3903 

Herndon, Va. 

142 

5 

- 33 

- 26 

4082 

Merritt Island, Fla. 

- 12 

3 

- 27 

- 23 

4280 

Vandenberg AFB, Cal. 

91 

3 

- 30 

- 32 

4050 

Pretoria, S. A. 

1604 

' 

6 

- 1 


4742 

Kauai, H.I. 

1157 

- 

9 

- 4 


7036 

Edinburg, Texas 

48 

3 

- 11 

- 12 

7037 

Columbia, Mo. 

249 

3 

- 20 

- 24 

7039 

Bermuda 

- 5 

3 

- 37 

- 36 

7040 

San Juan, P.R. 

9 

3 

- 41 

- 41 

7043 

Greenbelt, Md. 

27 

3 

• - 29 

-26 
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Table 2 (continued) 


Number 

Station 

Constraints 

N 

H 

a 

MPS7 

L 

'Vincent et al. , 1971] 

7045 

Denver, Col. 

1767 

3 

- 13 

- 13 

7072 

Jupiter, Fla. 

- 10 

- 

3 

- 26 

- 24 

7075 

Sudbury, Canada 

251 

3 

- 30 

- 31 

7070 

Kinsgton, Jamaica 

423 

3 

- 22 

- 23 

8009 

Delft, Holland 

72 

3 

48 

47 

8010 

Zimmerwald, Swiss. 

957 

3 

51 

54 

8011 

Malvern, England 

165 

5 

57 

52 

8015 

Haute Provence, Fr. 

702 

3 

59 


8019 

Nice, France 

432 

3 

45 

55 

8030 

Meudon, France 

214 

5 

47 


9001 

Organ Pass, N. M. 

1633 

3 

- 16 

- 18 

9002 

Pretoria, S.A. 

1564 

G 

0 


9004 

San Fernando, Spain 

81 

3 

45 

55 

9005 

Tokyo, Japan 

99 

6 

41 


9006 

Naini Tal, India 

1874 

8 

- 45 


9007 

Arcquipa, Peru 

2477 

9 

19 


9008 

Shiraz, Iran 

1588 

10 

- 20 

, 

9009 

Curacao, N. Antilles 

- 19 

5 

- 29 

- 26 

9010 

Jupiter, Fla. 

- 9 

3 

- 27 

- 24 

9011 

Villa Dolores, Arg. 

618 

8 

6 


9012 

Maui, Hawaii 

3036 

9 

19 


' 9021 

Mt. Hopkins, Ariz. 

2362 

3 

- 37 

- 22 

9028 

Addis Ababa, Ethiopia 

1911 

10 

52 


9029 

Natal, Brazil 

37 

10 

- 12 


9031 

Comodoro itivadavia.Arf 

;. 215 

15 

- 9 


9051 

Athens, Greece 

242 

5 

46 

54 

9091 

Dionysos, Greece 

454 

3 

78 

54 

9424 

Cold Lake, Canada 

684 

6 

- 30 


9425 

Edwards AFB, Cal. 

756 

3 

- 24 

- 28 

9426 

llarestua, Norway 

622 

3 

45 

46 

9427 

Johnston Island 

17 

10 

28 


9421 

Riga, Latvia 

32 

3 

22 

24 

9432 

Uzhgorod, USSR 

236 

3 

45 

47 


-Rft- 










one part in 750,000. The relative constraints used and their weights 
(l/cr 3 ) are all given in Table 3. 


2,44 The Adjustment 

The four sets of normal equations (See Section 2. 42), 
and the previously explained constraint equations were added together 
and a single solution was obtained for the combined systems. 

We decided to run three different adjustments to investigate the 
effects of the constraints we were using: MPS7 was ultimately 

chosen as the best adjustment, It contained all the constraints 
previously explained, inner adjustment plus height constraints. 

MPS8 included the height constraints but without inner adjustment. 
MPS9 was non with inner adjustment constraints but without holding 
the heights. 

After MPS7 was rim, we immediately computed the undulations 
(N) at selected stations and compared them with the values given in 
[Vincent, et al. , 1971]. This comparison is given in Table 2. There 
are some discrepancies, but generally the fit is good, indicating that 
despite the free adjustment, the height constraints had held (thus our 
origin is reasonably close to the center of mass). 

The results of the MPS7 adjustment are tabulated in Appendix 1. 
The number of degrees of freedom was 10586; the quadratic sum of 
all the residuals 12201; and the standard deviation of unit weight 1.07. 


2.45 Comparisons with other Solutions 

Table 4 summarizes the transformation parameters (systematic 
differences) between the MPS7 coordinates and those published in 
[Gaposchkin and Lambeck, 1970], and in [Marsh, etal., 1971], for 
the global network and for both the European and American nets. 
Two sets of parameters are listed. The first was obtained through 
the assumption that only translations exist between the sets of 
coordinates. In the second solution, the rotations were first computed 
through direction cosines independent of translations and scale factor. 
Subsequently the general seven-parameter transformation was carried 
out with the three rotation parameters constrained with their variance- 
covariances obtained in the direction cosine solution. Appendix 2 gives 
the general solution and variance-covariance and correlation coefficient 
matrices obtained in each case. 
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Relative Position Constraints 
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APPENDIX 1. 


(all coordinates in meters) 


ST AT I UN NUMBER - 
PREL. COORD. - 


CORRECTIONS — 
AD J. COORD. - 


102 1 • PLCSSCM FT . MO. .. ELLIPSOID 

ii 16002 iVifae ’ "’-4876023.0294 3943045.8335 


14.3103 -6.9144 -50.3740 

1118016 .7271 -4876329 .9438 " 3942995.4567 


VARIANCE-COVAKIANCF MATRIX OF THE STATION POSITION 


12.313516 1.496271 -C. 50642* 

1.49627 1; ... 9.574184 ..4.554761 

-0.5064.84 4.554761 11.640591 



STATION .NUMB E R 1022 


PORT MYERS 


_. ELLIPSOID 


PR EL. COORD. - 


CORRECTIONS - 
APJ . COORD* - 


. Y L 

8078 21. .32 34 -5o5 19 92 .6336 2533574.2911 


14.6979 -3.3427 -43.3542 

807646.0263 -5651996.4763 2633520.9369 


VARIANCE-COVARIANCE MATRIX .OF, THE STATION POSITION 


8.464229 0.425609 -0.775519 

. ... 0 i4 25.609.. 4.69176.4 2. 506412.. 

-0.775514 2.508412 8.352733 


STATION NUMBER ' 


PR EL. COORO. - 


CORRECTIONS - 
AOJ. COORD. - 


-23 57271. i-3 64 -4646320.9594 


13.2372 -10.9745 

-23572 57.3992 -46463 31 .9338} 


ELLI PSOID 

/. 

3668373. 0.24S 


-46.2868 .! 

3668326.7380 


VARIANCE-COVARIANCE MATRIX OF THE STATION POSITiION • 


6.934467 - II ,73.9 343 4.2451 14 

1.73334 3 I. _-.12 .95^55.4...:'. .1.1 ._ 1. 7.546.96 

4.245114 1.754896 12.444189 


STATION. NUMBER - . 

HH 

ST. JOHN'S IP, r 

f-P. . ELLIPSOID 

PREL. COORD. - 

• 

.. ... . 

2602652.1506 

-3419 24 0.9 2 2 <6 

4697706.79! 1 

CORRECTIONS - 

31. 1331 

15.08177 

-3 o??l 

ADJ. COORD. - 

26026 R 3. 253 7 

-39192 2 5.8409' 

4 697 vo9 • 8<> A 0 

VAR I ANC6-C0VAK 1 Ai\C 

E MATRIX OF TH 

E STATION PUSiTVlLN 


18 76. 78 
18 i 6.2265 132 
42 5.1 JJ3 3 6 

I 616 .22«5i>2 
2 2 •• 5 . V 2 8 02 7 
5 a 9. 0 30 2. 2 6 
-64-. _ ...... . . 

42 5 • 1 09 3 J6 

. . 569.83022a . 

19 9.52015.0 





STATION. NUMBER. - F A ! -^dA^KS ALASKA ELUPSUIO 

PREl . ’COUHO'. - 22^ ? 2 r3 LoP?. 7 " -1445669. 3178 ~575i 89^71 79 ' 


V^F’AI :• p.ArlKS -MASK A 


' CORRECTI ONS - . : 05 S 129.6797 -2 8 . 4 1 7 7 

/. AOJ.;tctpRi0^^- 5 ’-1445559. 63.81 '3751828^3003 

v, v a i ; /i:Ki <4 ^ ^ . j h ci ' ~ s't: a f ion "pgs i t.i on * 


STATION. MJ MB E K. 
PREL. COOkO. - 


j^k*. 1 0 ^;;/;.. . . • 

;'jv£ ■(■■. 

’'Wsiifil.l ’ 84?" 


r. GRAND PORKS MIN ..ELL I f : SC ID 

-424 20 49^7016 371876 5TV lTl 


AO J. COORD 


VA Rt ANCE.TCyV.AR IANCc/. MATRLX OF THE STATION POSIT I Of 1 


12.3V6699 

-3.443950 

-2.388826 

-3. 44395 0... .... 

11.786131 

6. 0340.13. 

-2-i 336*2 6 

6. 054 Cl- 3 

10.666619 


...£LL IPSO ID 


| C 0 R R E C TiT 0 N S - 
ADJ.^CbOROV: ^' / 

?:9935 
647*90. 998 3 

-3.9251 

-5177940.8080 

-44.9581 

3656736.8703 kt 

1 •■/• . Ofj-r 

':. V A R I.A N C E - C 0 V A R I AN.CE MATRIX/ OF THE 

STATION POSITION 

7 . '•’ 7 ■ ' 

•' 4 11 . '1,47>40(,' \ 

-0.1 2 7'4.6b 

-1. 203739 


40.127 : t : 4o : 

, .. .. 6 .986618 - . 

.; ...3, ‘‘40431. 

T : --: 

..'J;*l;. 5 2p373 : 9 ^ 

. * 3.440431 

1 0. 6. 5° 8 20 


r-CPRRECTIONS - / ,11. 3-26' -1.73*5 

' r \/: 2 3dli37:.05 1 3 -5372177.3368 

' VAR! ANCEHldy AR 1 AKCF I ATP I X OP THE ST ATI ON POS I T I ON 


—4 2 » 68 8 5 

1668571 . 6009 


-»..6l-«4>s7 2.290235 2. -350979 

,;3 ...; 2 i 2'»02:0 3 ... 8 . 5 30 703 .. . 2.24 4 506 


.:2V.'5 8y9 ; ? 5. ; 


3. 24° 50 1> 


15.227372 






ST AT 1 ON; 


J.LL.U’S.QUl. 


. .-'• ... Y ... . 

i-Kcl. CGflRft/. . , - -35.004., *23 4 -533 7967.2.543 3^*33514 i 22*^ 

- i ■•.■ ■■• • / . '/ 

CtJHKtC.Ti-Cn.S- >. ' •■> ■ 20*1273 . N -13.3o0b -57 .2096 

AO J . CCg f <L). - . -349H .*961 -532791^0152 3^9245 77 11 96 

VAR 1 ANCE-CCVAR 1 A-'.C7 NAT*.! X .OF THE. ST.i’.T I ON PC-S1 tYcn 


3 -> :: 
! f j . 4 t> 'f 3 1 3 


56,074 34c 
43.305331 


,43.50 58 31. 
8 5 .515 0 C 0 


STATION NUMBER - 


;CL.Cr..A PC SPRINGS. . 


.ELLIPSOID.. 


PREi. COORO. .- -1275219.15-.? -4793014.2 603 ~ 399^23 1 . 6200 


CORRECTIONS - -3.16 5.5 ' -11.9150 -46.6123 

ADJ. COORD. - -1275222 *3 197 -479 5026.1733 ~ *399423 3 . 00~1 7 ' ’ ~ 

VARIANCE-COVARIANCE M AT K I X OF T HE STATION. POSITION^ 


.... 130.522 156 5.439216 1.453025 

...; .. 5.439iT6 :... ..._ ..39.43172,0 ...18.06 7636.. 

1.453025 18.067686 23.139421 


ST AT I ON ..NUN OCR.. - 



PREL 7 C0CkB 7" - lVi ? 12 T. 14 53 " " ’ -4 4635747V67 7 42621: 


CORK EC T1 C.'iS - b . 3 5 R 5 

-5.3602 

-50.9707 

AO J . CCUF-O. - 151 3129.501 S 

->463530.0620 

42330b 5 . 2927 

VARIANCE-COVARIANCE MATRIX OF THE 

STATIUN POSITION 


14.199741 

2 .9 5604 2 

0.642642 

... 2 • 95604 2 

, 19.573991 

3 . 80 l -7o5 . 

0 • 64i6-t2 

3.808765 

14. 255383 


ST A I ION NUMBER - 

PREL. COORD. - 
CORRECTIONS - 
ADJ. COORD. - 


SEMM3S A i.A * A u 


.t.LL.IES.U.1.0. 


. X ... Y. ;..... _i .... 

1657235.3703 -644 1576.7655 224511*. 916 = 

!o. 4-54 5 -1,1263 -47.196-8 

*47252.3 25C -5-.31977.K958 •” 3 7.4506 7.7200 


VAR IANCE-COV/iRl ANC.- fATKI X T hi STATION POSITION 


1 8 . 0 1 -3 3 :> t 
-0.022277 
0.716 3-'. 


-0.02727? 

. 8 . 1 v ? 31 4 
4 .274106 


0.71 n049 
.. 4 . 7 7 i 06 
14.633314 






S T A T I ON N U f* y 2 A - 

''.3404 

SWAN ISLAND 

. ELLIPSOID 

PREL. COCAO. 

' X 

44a 47( .3473 

. _ Y 

-605 39 59.6531 

189570°. 1149 

CORRECTIONS - 

r: • ■> 5 

22.5346 

-55.5610 

ADJ. CCORD. - 

r>424T0._2358 

-6053537.3184 

) 39571 3.5539 

VARI ANCE -CO VAR 

ICE MATRIX ‘Jr THE 

STATION °C S I T ! CN 



26.22159 8 
. -1.935E13 

-1.933513 
2 *.,00356.6 

5 • 3004 C 9 
.. 0.7i!tt6 3? ... .. 

■s- 

5. 500409 

0.736632 

27. 69 i 902 

STATION NUMBER -. 

.. 14.05 

DP AMD TURK 8.1. 

ELLIPSOID. 


X 

... ,y._ 

7 

PREL. COORD. - 

191° 4 74. 3507 

-562 1096 .3563 

2315643.4691 

CORRECTIONS - 

h . 5264 

-9.5606 

-29.2670 


ADJ. CCORD. - 1919479.3772 -5621107 .9139 2316814.2021 

VARIANCE-COVARIANCE MATRIX OF THE STATION POSITION 



- 10.192175 
2.84 76 )3 . ... 

2.847803 

10.2.1.1.54,4 

3.290664 

3 .5.3 1 1.7.1 _ .. 


3.29066-4 

3.561171 

14.633133 

" 



i i ii mum 

STATION NUMBER - 

3406 

C.m LC. AO N. A NT III 

JL1.- _r.U..!P SO 11). 


. . .-. X 

,• Y 


pr-FL., CC’CRD. - 

2251788.8701 

-5316921 .9460 

1327276.8528 

CU^ CC T I Li M $ ■ — 

5 . :> 5 7 7 

-0.7047 

-53. 7261 

ADJ* CC'JfcD • - 

226179 o«4? 64 

-53165 22 .6527 

1327223. 126 P. 

VAR I ANCE-Cl’V.V- 1 ANC 

- >*• ATRIX i : f Tri 

= STATION POSITION 



11.255736 
... 2 . 0 2 i 9 i 1 

-2. 4 8 SO 9 7 

2. .02391 i 

. ... .6 .224379 

.0.212203 

-2.466057 

... 0.21.22 0.3 

16.691533 

STATION NUMBER - 

■HHH 


„ 'ELLIPSOID 


X 

y _ 

. . 2. 

PRSL. COCRO. - 

2770863.5463 

-5513555.9 595 

1131225.6195 

CORRECTIONS - 

-2.3249 

18.3334 

-55.103;-: 

ADJ. CCCKO. - 

2779881.31 *9 

~-oE13537.57c.l 

11.61 174.516? 

VAR! ANCfc-COVAR I ANC 

E MATRIX OP TH 

E STATION PCS IT ION 



37.511 4 i 5 
3. °69ri 0 
- 2 •• . 4 1 K 5 1 0 

8 .9 35 610 
. 2 0.4 55 6 9;. 
67 _ .6 5440 c 

-24.4 i- 510 
- 14. E 5 '-4-J3 
67.3-16878. 




STATUS NUMBER - 364S 



Pr EL . CJJRP. - 
CORRECT I liNS - 
ADJ. CHORD. - 


832544.6662 -5349: 

1 2. 2 7oc -1. 

3 ‘ 2537 . 9^-30 - 53 -°: 


- !:LLI P.SQI D 

54 3360660. >}i’f 9 



VAr.I AriC t-CCiVA:'. ! ANCE V.TIUX Of THfc STATION POSITION 



ST AT I ON NUMBER. - 
PR EL. COORD. - 


CORRECTIONS - 
ADJ. CCOP.L) . - 


— 1 . 4 H 6 A 3 1 
,_fc. 722818 . 
2:0.223213 


1 ABERDEEN' .".ARYL AND . . ELLIPSOID 

li 86768. 5113* -4 785190'”. 10 73 ” *4032960.5235 


12.7733 -10.1142 -46.2572 

1186731 .2646 -4765200.2215 4032914.2533 


VARI ANCE-COVAR ! ANCF MATE IX OP THE STATION POSITION 


- 13.786196 2.626266 -0.030759 

.2.6 36 26 6,. 1 2. ,.8.6.920.6 5. 30.6.2 4.-< 

-0.070759 5.306 246 12.63.4535 


SEAT ION N'JMcHk 
PP-fL. COG 2D. - 


CORRECTIONS - 
ADJ. COORD. - 


- 3.3 M HCHgSTi AU, P.LQP I DA _Ei L I £5.0.11) 

961749.7864- -5679161 . 02 94 27 2^9 6 0. *9 76? 


12.57*7 -3.1751 -45.332o 

961764.4651 -5o791 64.2645 2 7299 15. 643 7 ~ 


VAR I A N'CE-C UV AS I AMCC MATRIX .JF THE STATION POSITION 


12.21 7638 
0.640301 
-1.5 6 69 ‘JO 


0,630201 -1.55=960 

. 6 . 1 3.1.9.25 .3.4.43 

2 . l*»li «4S 11.199167 


STATION NUMBER.- j 

PREL. COCRD. - -12: 








Jtm 4 ■ ? i — 

ASO.I 

HERRON p v 100 

U .. ELLIFSUU). 


X . 

. . . V.. 


PPEL. COONC. - 

10S8973.44-S2 

-4343050.3590 

3° 91 805 . 25 9 C 

CuKKEC t IONS - 

2.2447 

. 38.3767 

1 . 9134 

AUJ. CCOr.O* - 

1 0 H 9 0 7 5 . a 9 7 V 

-*t Ot 30 1 1 . 9 S2 3 

3991811. 1714 

VAkI AiSCc-CwVAs I A.'JCt 

MATRIX IP THL 

STATION POSIT I CN 



•152;.*5?4<.2 
. 5 . *"• 09G-. *• 
-TO.. 15 34 6 1 

5.903044 

‘70.053888 

4-3.60X731 

-30.158461 

...45 .6047.21 . 

62.065319 

STATION NUMfcuk - 

4oso 

PRETCR! A 

__ ._£ O.L.LP.SiJ I. IX 


• X 

V 

.... . i 

P R E.l COORD. - 

•SOSlfcOS. 7772 

2726595.1416 

-2774193. 7107 

CORRECTIONS - 

42.5116 

16.5502 

97.8050 

ADJ. C-C-ORD. - 

50 5,ifc62. ?P£9 

2 726.511.7219 

-2774095. 9007 

VAR1 AVCE-CDVAR I AN CD 

MATRIX ih Thu 

STATION POSITION 


• - •• ’ • 

1 1 4 . 7 0 7 7 '> 7 
174. toO 'rl c 
343. /.Oi l 40 

174,660583 

. . .6.5.7 *.4ii, 6.42 

1049.241343 

342.2011 60 
1049. .24.1 34.3. 
1 7 5 2 • 7 4 >- *, to 

STATION »\ U . v > 8 c - 

406 L 

ETRANT 

ullips'cijo.. 


. /. X 


J 

P P EL V COCKU. - 

2 88; 594. 3 640 

-5372.524 .5483 

1660020.2732 

ccpkCct IONS - 

-3.2296 

-12.2235 

3fc. 9964 

AUJ. CCJKU • ** 

2331591.13-4 

-5372536.7716 

1866057.3717 

VAR ! ANCE-CUVAR I-ANCc 

MATRIX OF 

STATION POSITION 



6 . £• 1 42 3 4 
2.2 9UJ8 3 

2.290063 
.ii.5 3l6.73 

2. 551892 
3.249.712 .. . 


2.551 592 

2.249712 

15.222010 

STATION NUMBER - 

nm 

Tkf.RT 

.■XLLIP.S.OI.O 

PR-fcL « COORD. •- 

X 

19 204 J l. >J65 A 

Y ... . ... 

-56153 96.5502 

Z 

2219115.5044 

CORRECTIONS - 

5.9831 

-3°. 0562 

51.4104 

AO J. CCORO. - 

19 20.4 0 7.8440 

-561 >37. 6064 

2314166.9142 

V/.RI ANCc-CUVARJ AMCc 

MAU I a OF THE 

'STATION POSITION 




10. l'-io+c 
. ..2*64807 3 
3.300149 

2.346072 

,10.219 572 

3.580917 

2. 300 i 69 

3. 5R0'H7 

1 4 .55 7j b I 


*' 69 “~ -- ------ 





ST AT I UN iVJNBcR - 
PRBL. CUCr’D.- 

corrections - 

Ab J • C 00 is 0 . — 


910571. 7663 

— i 0 . z 6 i> v 

910560. 510? 


YAklANCE-COVAP •■••iTF. !x OF T: r 


. 0. .98sy:)I. 
-0 . BA '■•06 2 


-062 6109. SI f. 2 * ? 0 1 7 972 .26 2 ? 

- 10.36-9 2 3. 9777 

-053911V. 6631 3017905 ,7'*05 

STATION POSITION’ 


•986981 — 0 . 6 A 6 5 6 2 

2 . 'Vrb.l.o.r?. ..2 • 3.205.5o 

2. B 305 06 7 • 7 1> 2 6 A 9 


STATION NUMBER - 
PREL. COORD. - 


CORRECTIONS - 


• ... X 

•207385?. 93 


-15.2215 


-452 12 C4 . 57 A 


-6.9001 


-a. 5 2 A 2 


ADJ. COORD. - -2671675.1 595 -A 52 12 i 1 . 47*9 2607001.1501 

VARIANCE-COVARIANCE MATRIX OF THE STATION POSITION 




1 A . 7 5 A 6 5 3 

0.959315 



... 0.95931 5 _ . 

1.6 .04.0498. __i. 



5.7-849 5 

15.904103 2 



station number .- 
PFFL. COCPO. - 


■CORRECTIONS - 
ADJ. COORD. - 


4740 NR ER 34 _£LU. P.SUIJ1 . 

2306892.0318 -467 4293.01401 330303*-. 91 34 


-10.4465 


-15.269? 


2308861.5853 -4674306.282? 


30 . 26A5 
12931 i5. 17 7 b 


VARIANCE-COVARIANCE MATRIX op Tile STATION POSITlCN 



32*40 

-1 .34 391 9 

— 4.421 1 4 4 

.*> § i s 

9 .9 3153.6 

3 . c 4,6.7 43 

11/* 4 

’5.64674 2 

15.94000 9 


STATION NUMBER - 


mm 


..ELLIPSOID 


.... Y /. 

205 5 -9 • AO 5“"' 2387495.291} 6 


PkEl. COORD. - -53 


CORRECTIONS - 5.55<i5 18.830*, 11.6629 

ADJ. COORD. - - 55 / »39?2. - 72 5 —20.64530.52 53 23 87 506.9616 

VAR I ANCE-CUVAR I AN'Cc MAT*. I X Jf THE STATION PUSITilON 


1 6 . 2 j r 
8 . vt>2 • i 2 
-2.07135- 


h . ° A 2 9 n .? 
.. 23. 169 >37 
-1 .645237 


- 3. 0? 13=«. 

- 1 . «'» 5 ? 6 7 
3 2 • 8 1 * 5 2 7 







STATION .NUMBER . -... 
RREL. COORD.' - 


CIILNSURC TjjXAI 


cLL.I PSOi U. 


-6P.b519.5T2 7 -5 657465.9567 2816385 .7705 


CORRECTIONS - 

23.2528 

-7.7280 

-43. 7482 

ADJ. CCURD. - 

-•J 26496. 2290 

-5657473.6647 

2516342.0224 

VARI ANCE-CUVAR 1 ANCE 

MATRIX OR. THE 

STATION POSITICN 



14 • 34 »i4 0 

-0.792797 

1.53-558 


. -0:. 7.92 797 _ 

..2.8.7 2.3 8 . .... 

.. 2.. .1.7.3 50.0. 

. - 

I". 52 255 3 

2.173500 

12. 162002 

STATION NUMBER _r 

7037 

COLUMBIA MCI. 

.....ELLIPSOID. 


X 

Y 

, 

PR EL,. COORD. - 

- 19131 7 . 0 or. 4 

-4567282.1392 

358332 5 .7103 

CORRECTIONS - 

lb. 1420 

-14.2615 

-44.1249 

ADJ. COORD. - 

-19129.8. 9264 

-4967296.4507 

3983251 . 5854 

VARIANCE-COVARIANCE 

MATRIX OF THE 

STATION POSITION 


* 

11.105130 

-1 .437472 

-1.437143 


-V.437'**7? 

6.79 v r» 9 0 

3. 6524 3« 


-1.437143 

3.652438 

9. 149429 


S T A T. I GN_. NU tf 6 6 
. P R E L •” COGR D 


CORRECTION'S - 
ADJ. COORD. - 


2306200.5140 -4873602.7626 2394650.8205 


6.3249 -3.9950 -39.5362 

2306207. 3589 -4873607.7576 3394591.2643 


VARIANCE-COVARIANCE MATRIX JF THE STATION POSITION 


9.304137 — 1.34 7172 — 4.419626 

.. ....... — . ..’ -1.34 717? 3 .92,9902 ,84 7 7 07 ,. 

-4.4iq.s3fr 2.347707 15.934007 


S T.A T I ON NO M HER . - . _• 7Q~p 


. . SAN JIJA.N P . i! . 


... .ELL! PSOIO 


PP-EL. COORD. - 2*65036.6657 - 552**931. 4949 1965565. 1070 


CORRECTIONS - 8.4647 -9.2015 -33.5274 
AOJ. C'luRJ. - : ; 2466045. 1304 - 5534540.6964 ] 5o5S* 6. is 3s 
VARI ANCE-CCVAk I AMCE MATRIX OF THE STATION POSITION 


16.731550 
4 • o ) ’«? V t 
- s. 2 J 3 5 i 6 


4.01574s 
>• . 096200 . 

U . I n v 6 5 


-6. 33 3528 

0. 1 ••..? n r.8 

17.420575 





ST AT I ON. .NUMBER _ r 


70^3 , GREENBELT - v A k’ Y L A N D EL LI. PS 0.1 0 

. X Y 2 

1 130691; 30-35 "4831330.1355 3994?! 1 . 7359 


PREL . COORD. - 1 130691; 3035 -4831330.1355 3994?! 1.7359 


CGP.k cCT 1GNS - 11.3710 -5.3198 -46.21 72 

AO J. COOP.O. - 1130702 • 5795 -4831336.4753 359M65.5687 

VAR 1 ANCE^COV AA 1 ANCE -LA Tk IX GIF THE STATION POSITION 


13.00^326 • 0.72b 5 54 -0.S7S7C2 

0 . 7 2 -j 9 34 9_. 2.3.1 2 26 2,3155.93 

-0.975702 2. 919598 5.754oS3 


STATION .NUMBER.-. 


. DEI-iVE K COLOR ADO. E 


PREL. COORD. - -1240451. 208 3 -4 760227.8963 ^04905^.2245 


CORRECTIONS - 12.9239 -7.1289 -48.8968 
ADJ. COORD. - - 124 04 78. 2.84 4 --*7602 3 5.02 52 4C49005.3276 
VAkJ ANCE-CCiVAP I ANCE MATRIX UP THE STATION POSITION 


19.376439 • -3.993436 -0.808182 

-. 3 . .9 9 3 4 3.6 1 0 . <>.?. 3 798 ;■ 4 ... 6 . 3.6 1 . .9 3 

, -0.808182 4.638193 11.860496 



STATION .NUMBER - .. 

7075 

SU'D 8 OR Y CANADA 

_ ....ELLIPSOID 

PREL. COOKOl r ' 

.. X 

6925-92 .3702 

Y . _ . 

-4347069.9512 

2 

4600546.1501 

CORRECTIONS - 

?0 . 1 09 2 

-8.5565 

—4 1 . 4 j 8 4 

ADJ. CUORO. - 

4 r 2 61 2.9 795 

— ,347076 .509 7 

4 600 504 .6-517 

VAR I ANCE-GOVAR 1 ANCt 

•MAT 9 I X OR 7 

STATION POSITION 



17.116149 
1. 828 309 
-0.446265 

. . . ■ -•* ” 

1 .5 23 30.9 

18.262697., 

9.04 1 JOti 

-72- 

— 0 . 4 4 ■* 2 6 -J 
. 9.04 15 3-1 
15.230182 












STATION .NUMtfSR.-. --.- -3015 HAUTE PROVENCE. F !•". ELLIPSOID 


_ _ t . y _ _ 2 

PREL. GOCRO. - 4578327. 1964 43 79 63. 9 9 36 4 403179. *0029’ 


CORK 

fcCTICNS - 


-0.3565 

-5.9354 


'47.6372, 


AC J. 

COORD. - 


4576327. 6399 

457960.0432 

4 0 3 2 2 6 . 6 9 0 0 


VAF.I 

ANCE-CCVAR 

I ANCE 

MATRIX UP THE 

STATION POSITION 







23.86565c 

' 1 .970908 

-16. 

39-:-.?5 





.... 1. 970905 

_55. 757516 

■* 6 * 

6 1 .2 0 ° a 





-16. 36r>335 

-6.612006 

24. 

45 5'* 00 


ST AT 

I ON... NUMB ER 


8019 

NICE, ...FRANCE. . 

... .... 

HI! P.S:: 

ID 




X 

Y... 


..._i 


PREL 

. COORD. — 


6579465. 99/9 

566599.0005 

4 3 36 

406.0076 


I 

o 

o 

/V 

EC T IONS - 


2 . 3669 

-8.1364 


44 .4101 


ADJ. 

COORD. - 


4579466. 364 7 

5 6 6590 . 564 1 

4 366 

>452.4177 


VAR I 

ANCE -CO VAR 

I ANCE 

•MATRIX Of- THE 

STATION POSITION 







21. 93749 A 

0.151565 

-17. 

395617 





9.151965 

49 .2353 76 

..-5, 

219443 


| 



-17.395617 

-5.219443 

• 

21 . 

593170 


sr«r 

i ON NUKStrR 

- 

SO 3 O' . 

« ROD ON FRANCE 

M. 

xLLI.PS.D 

J.D.. 

I. 



.. X - 

_ Y 


2L. .... _... 


PR El 

. LGuRD. - 


A 20 5 3. 8:651 

163731.9484 

4 776 

53 1 . 84 0 ?. 


CORK 

EC T IONS - ■ 


• 4 50.262 2 

-26.1516 


-9.0951 


ADJ. 

C CORO. 4 ; 


420563 3. 9 82 9 

•. 1 . 

163705.7963 

4 776 

5*7 2. "*74 51 " 


> 

7? 

M 

ANCE -’CO VaR 

I ANCE 

MATRIX OF THE 

STATION POSITION 







143.494990 

33.360690 -103. 

562842 



; • • , 


. 3 J.i!b06 rj 0 

179 .540 <6 5 

-7.2. 

2.94619 




‘ ** 

1 0 3 . 3 6 i 6 4 c 

-72.29661 v 

124. 

032169 


ST AT 

ION MJM4ER 

- 

9001 

ORGAN PASS, In.M; 

- 

ELL IPSO 

! 

I D 

PFEL 

. COOkC. - 

- 

.. . . X 

1535756.9990 

- 5 1 6 c 9 9 5 . 9 9 5 4 ~ 

3401 

. . Z._. 

04 2 . 00 3 3 

--- 

CCPR 

ECTluNS - 


1.157b 

-14.3514 


2 1 . Oo77 


ADJ. 

COORD. - 

- 

1535755. 1114 

— 53 a 7 010 . i?46>« 

’3401 

063. 0710 


VARI 

ANCE-CCVAR 

I ANCc 

MATRIX OF THE 

STATION POSITiON 



1 


21.67 343/. -5.362339 
-6.8623S9 6.45690* 
-1.9 120?!, 3.044096 


-1.632074 
3. 06 •'•->96 
10. 7o-*7 c7 





ST AT !JN NO *8 £9 - 


PR E TORI 


PP CL. CHORD. - 


COR* ECTIlNS - 
AD J. COORD. - 


. . .... X Y 

50-3612 5.0005 2716511. 00:7 -2775784 .0066 


'7.2.306 6 • 2 d 6 6 66.0435 

30 56 132.2110 271651 Y.Z* 1 2 -7773 & 97 .9630 


VARIANCE— COVAR-IAhCE rtATKlX OF Tri- STATION PCS IT! CM 


1 1*.. 68 2 3 -? 4 
.17*3.62 60.' 7 

5 A 3,. 1 7Uo oi 


STATION ..NMJKijCK jr.. 

5004 SAN Pc' 


X 

PREL. COORD. - 

51 05 5 67 .9960 -55522 

I'- - .. 


1 74 .o 260 37 343.17 06 6 3 

657.421226 _ 1049. 1 7.0 2 72 

1049. 1 702-72 1 7 5 2 . 7 <’* 1 1 53 


SAN FERNANDO, SPA! EL L.I.P.S.O.i.!). 


CORRECT IONS - -4.9963 -21.3000 **2.4355 
AOJ. COORD. - 510,5532.4992 -555249.2010 3769709. -*926 
VARlANCE-CCVARiANCT MATRIX" oT TH£ STATION "POS iYlciT" 


- 13.135126 -5.541513 -10.721963 

..... -5. .5.4 151 3 114.2.02403. . 11. 0 3. « 4.7. 5 

-10.721963 ' 11.038475 17.726254 


STATION. NUM6-SR 
PR El- COORD. - 


"CORRECT I DNS 
AOJ. -COORD. - 


TOKYO, JAPAN 


E.LLI.P.S.C’.I.D. 


-3946692.9973 3366253.9960 3698852.0033 


1 6 . ' ) 4 6 C 


32.2762 


1 fi. 73 2 E 


43946676. ->4«2 23663 31 .2762 3693250.7366 


VARIANCE— CoVARI ANCii MATRIX >iF THE STATION POSITION 


:.!«>! 0.-62541 2 1135 .7 67 243 370.520660 

3 9.7672 *3 1 3 i 5 .3.3 3.91 0 190. I 4 3.97 o... 

7U.52u--.60 190.14357c 242.705340 




CORRECT JONS - 


-42.9171 


55 7 1 1 02 .996 2 


10.3435 


33.0793 


AOJ. COORD. - i 1013160.0815 54711 13 .339? 31 09661 .0914 

VARiANCc-CrVAR IANCr MATP IX Ur THE STAT ioN POS i TICN 


3 52.->1''7 )•; -121.490997 

. ...... !- 12 1.5409 >? 99. 3 5n 59.9 

5°. 6 7 75 >*-. _75_ —33.65132s 


5 5 . 7; 7 7 ft 5 i 
- 17.' -.-i 3 2* 
76. 26.2703 







!J.ILRS&IU. 


X .... . ...Y ...L . 

P'f El . CO .nppv. 194..: 7 74.9 ; t98 -5 zO-iO 11 .0013 - 1790923. 0040 

CO^tC-TICNS' -' •-•15.7,647 -13.4414 45.659c 

Ai)J. CGORp. - 1943759.3311 - 5 80 4 0 94 . 3 ft 23 -1796367. 14*t5 

VAF 1 ANCE'— Cl : V-A!VI ANCE .M AT W-l-x <JF THE STATION POSITION . 


.701490 
. 24 1<- >7 

.2210-4 


-0.241407 
. 1 7 . R 4.9.9. o 4 
-13.1 13526 


ST AT I ON NUNbER . 90 04 S H I X A 

PKcL. COOKD.”- 3376-892. 6997 ' 4.403975. 99-85 3136250.0045 


CORRECTIONS - -13.3703 

ADJ. COORD. - 3376679.6294 

2.6709 

4403973.8692 

43.4706 
21 36292 . 47 5 1 

VARIANCE -Vo V A fi I ANCE MATRIX OF, THE 

STATION POSITION 

76 . 802336 

-36 . 4 36143 

25.409600 

•7 ■•-■'•'• •••'. - 3 F. 4 36 143 

_76..A 4829.0.. „ 

-19.501727 . 

. :'■>:•• 25.409.400 ■ 

-19.501727 

64.66240 9 


ST AT f :ji\ N'JKrVf: 4 - 
PPCL. COOK J.: - ", 


COFk EC TICNS- ' 
AOJ. COORD.' : 


. .. x .94 . 

.2 25 16 25. 9 99 i -5»lo9l9.9982 1327160. 00? 5 


-7075-707 -6 .6 544 25.1242 

225; cOa.4'254 -5616925.5527 1327195. 1263 


'VARI ANCE-CWaR:! ANCC HAT R 1 X Of THE STATION POSITION 


11.2,5573:4 2.023912 -2.438057 

. 2.0? 3912 . . fc .224 379 0.21220? 

-2.44405 7 0 . 2 i 2 202 16.691522 


ST AT ION NUMBER. ,901 


PRCL. COOKC. - 


corrections - 

AOJ. COORD . - 


1UPITE2, FLA. 


ELLIPSOID 


. .. f . : X : Y j _:... . Z... : .. .. 

776290 .9 99 5 -5601396.0027 2S30240.C056 


-19.7314 -10.1094 24,250? 

.9 70.2 71 .2661 -56 0 140.1 .1122 2800264.2559 


VAK I ANCE’-CCVAa ! ANCe MATRIX OF THE STATION POSITION 


7. 3501 1 1 
0.93611 • 
-0. E476I 5 


O . ■) :>n l i 1 - 0 . '"-761 ■ J 

3 .841 332 .. . 2.32 7 7 >*■ 9 

2 . V 2 7 7 8 9 7.756 5 01 






901 1 

VI It./ DUl.UR E S t 

AF.G ELLIPSOID 

X 

. .V . 

. .... ... . L. . 

2230588.9964 

-4914572.9996 

-3355426.0070 

-10.9111 

-2.5516 

41.6426 

2280573. 96 73 

-49 1 •* 575 . 55 14 

- 3355334.3644 


PFeL. COORD. - 22>S0588'.9~..j4 -4914572.9996 

"c C RS*fcCT ION S~ — I r_ - -’iu7 nil ' -2.551tT~ 

ACJ. CCORO. - 22 SO 573 .9672 -491 4 57b . 5514 - 

VAR I ANCE-CCVAP i/.NCE f-1A7?! X OF THE STATION POSITION 


It*. .531 5 A 0 5.167772 -12.638693 

5 . 16777 C- ...53 .070,05,5 ,_„_-7,C. fr6 0.6.82___._ 

•viS. SfiS6V3 -70.680632 133.616794 


STATION NUMBER. 

4 x 

PRlL. COORD. - -5»66052. 7998 -2404282.000 


CORRECTIONS - -32.8294 17.0400 6.9332 
ADJ. COORD. - -5466065.3 292 -2404264.9602 2242179.9569 
VARIANCE-COVARIANCE MATRIX Or THE STATION POSITION 


10.657097 

5.381140 

-2.308945 

5.831140 

2.1 .4134 78 

. — 4 . 89.43 4 3. 

-2.808)46 

-4.396343 

• 

30.-730210 


-64. 91 


3 0777 0 


27.6650 


LLUPS.'J 


i 9 2 6. 00 4 / 


46.430 


ADJ. CCGRD. - -1936646.3650 -5077676.3298 3331962.434 

VARIANCE-COVARIANCE MATRIX OF THE STATION - POSITION 


524.632305 - 346 , 13 5 962 -213.696170 

-345.1359-2 .. . 26 1 .1 16.45.4-. 1.1.3..431.425 

-2 13.696170 175.431425 1 4 i . 9 2 9 3 6 1 


902 S MiMS ' ARAHA 1 , FI Tm ! 0 ... liLLIP 

45 03 750.0202 39o5?00. 9969 363372.0027 

- 1 r- 8 5 -4.5775 4 9 . 64 4 

4° 03753. 1 Ocj 7 3965196.4214 96392 1.6471 


VAMA\CE-r.i VAR lANCt. V.Tft ! X ‘Jr fit STATIUN PCSITION 


• 1 >v ' 0 * 
-.r. 'JOOO-: 5 

• F . r- 7 6 5-- ■; 


-w .030003 - n 7 3 

75.1 7101. r 3 3 . 8 4 : 6 6 

-3.5 .846265 190.710875 


l/> 










ulLI PSOIL' 


station isuK^fK - 


MiilALi 


PReL. COG'S D. - . )O il -3653556.0041 -6.5432 5.0004 


CL-KR F C T I O..M S - 
4DJ. COORD. - 


9.5709 


•26.6133 


3 . 2 5 74 


51 86470. 5 7** 1 -3653fe72 .6174 -654316.6430 


VAR.I AnCE-CIIVARI ANC? MATRIX Of- THir STATION POSITION. 


I'-O. 412546 
“i, 07 8 n 4 -• 


ST AT I ON KU^bsK' - 

PPEL . COCRiO. - 
CGRP EC f I Uw.S — . 


C C MS » 5 RO RJVACSVIA. ....... SL-.U.PS0UJ 


X Y. l ...^ 

16 . 93 * 02 . 9 . 99.0 - 4112327.9956 - 45566 - 9.0032 


•15.6 rO v 


-10.1655 


-3.5227 


AC J. COORD. - 1653737.42 i,C. -4122336.1616 -4556657.5259 

YAP I ANCt-COVAR lANCc .•! A T R i X Of- ITif. ST A T • ON .POSITION 


67.9956 3 


■STATION NUtXOcft 

' M5] I.:.... 

ATM, Nl.GRtlXF 

i 

_ £L.L 1 PSOI.ii ~ 


. x •■:'■■ .. 


; j 

PR EL. COORD. - 

4-606 8^7.341 7 

2085885.1016 

3903 534. 303 3 £ 

CORRECT IONS - 

. . 13.8M 1 

2-2.0139 

r 

5- 

O' 

ADJ, COORD.. - . 

46 06573. 2028 

2029713.1205 

3903 605.3 970 

% V A R 1 A K C E - C O V Aiv ! 

ANCF MATRIX OF. ,THs 

STATION POSITION 

il 

•- " •* ' 

4 7 .76-31.1 

-75.7 76157 

0. 1222 17 | 

; 

-76. 7 76107 

. .3.32., 3 3.5.4 5.1 

..“.5.2 ,.3 4 0A21 | 


0.12221 7 

— .5 2 .340421 

47.224241 | 


57 A T I l \ IsuMb r ' - 


P 3 EL • COORD; - 4 5 9615-7 ,001 / 2D;vA2'-. 5953 3912650. 00 -SO 

CGPRLCT1CNS - -10.6! 51 4.7407 40,7523 

AIJJ. COORD. — '45 »5146.:H?.6 203942-7.7402 351 2 6. VO. 7605 

: VA.nA-vCE-d0VAR l*i«CC ’•< A r* ! A Or T’-F STATION POSITION • 

2 3.7516 6 " -7. <Vj So*.* 


*.* [ : j v S OF , QE : •■• C •_ 

Y.. 

20; 5 3 


..i.LL.1 PS DI O 


'3.7.166 
-2 <.«* 4:j«: 4 
— 7 . •« r .3 V ■: 


- t A.>4wal«. 

105 .4 49 001 


A . • • 


» *♦ 



.99 f’ 


STATION NUMBER - 


2±2h~ 


CO LO 


Pr FL . ccomt. - 

CORRECTIONS - 


-126^837.9987 -3<-o6?.~ 3 

4.1202 


-11.6673 


3 3 6 


alj. c. cord. - - 

1264633. -3785 

—34 00395.6655 

5185470. 7415 

V A li. I A ft C. E - C 0 V A.9 1 A Si C E 

MATP! X 3r The 

STATION POSITION 



27.601571 

0.7.50156 

6.7 50 136 
42. 163804 

1.095511 

14.351202 

... '. ■ 

1.095511 

14.351202 

'i5«66-. 1 72 

STATION NUMBER - 

9425 

ROSAMUND. 0 ■■ L ... 

.. ELL.IP.SaiO..' 



y 

2_ 

75 RE L i - CUCNp. - 

2456010. 997 3 

-4624-420.9955 

3635035.0032 

CORK'tC T I ON S - 

-4.64/7 

-10.3622 

12.1509 

AOJ. cccrd. - 

24 300 IS. 844 5 

-46244 31 . 3577 

3635047. 1:541 

V A R I ANC 5 - C 0 V A* I AN C ~ 

MATRIX Or THE 

STATION POSITION 



-16.600406 

0.0078-32 

0.007884 

1 5.3.530.03 

6.0?67 9 3 

1.5... 2 187 3,2 


6. 066793 


13.216732 


28.766095 


STATION HUMBER 


9417. 


HAR.ES TO. \ . N’ORMY 


FLU PS 0.1 0 


...-.- x 

PR EL. CO On 0. - 3121-279.9 959 


y. 

6 5264?. 9993 


"5 512 70 1 . 00 14 


CORRECTIONS - . -13.1041 -.16.7612 64.2684 
AO J,. C 008.9. - 3121 2 66 . 8616 592626.2361 5512 735.2396 
VAR 1 ANCF-COVAX 1 AfiCr *ATF j A ir The STATION POSITION 


96 .£43751 
.2 1 .5577 09 
-55.23 >439 


21.567705 -55. 2 3343 9 

9 2.853 0.7 6. - 2 0.24 1 66 7 

-2 0.241 so 7 44.135516 


STATION NUMBER - 


942 7 


S ; 4Tn,\' ISL., PA C 


9 R c L • CCCRO. - 
Cf!N r < £C TIj.NS - 
AOJ. CUORi) . - 


a - Y 

■60U74C2.OO0S -1 11 1K55.99E4 


. c L.L I P S 0 1 u 


1 o2*.'73u. 0040 


-53.049€ 76.4625 -12.4444 

- 6007- 3 E. 0506 -111 17 83 .5359 162571 7.5590 

V.Aft Ifi HC S'-C C VAR 1 ANC c MA T'H i X lr TiSt ST /• T iCN POSIT 1C.\ 


,10 I .J4.V>(> 7 
-1*5. 5-6 / 7 
1-4 .40 7\6>i 


-10 3 .5 4*; 72 7 
4 6 7 .723 4 -r . 


- 79 - 


— } 7 1 .50 *'v*t 0 


5 9 ... 0 7 8 o 0 
_ - 1 71 . 50 c. 0*0 
1 03.414781 




S T AT I ON NOMcix* ' - . . $43 1 

.... .. A 

?■? cL. GtlOKD.. - 3i ?3901 ; .03i5 


-.CORRECTIONS - 

• -l .529 i 

-3.0555 

63. 5594 

AOJ. C COP D . - 

31 R3eW'. 5033 

1 *t 'i !»r^4 • Co U 

5322335.5334 

VAK! ANCfc-CCVAN I 

AN CL •« A TP IX Or THE 

S T A T ! L' i\ POSITION 



• 153.746692 

-65.0 70518 

-59,5777 0 4 


.' ... -o.5.07o;91S? 

73.9 34434 

22 . / *.< 5.0 51 


-95.9777o4 

22 .42505! 

64 . 733.047 

ST AT! OK NUM&ER 

943 2 

UZOHORCiO t US S P 

_ tlU PS.CID. 


. .. - -V. x 

Y 

7 

PR£t; CC’Of'D;. - 

. 3907420.994? 

1602397 .0020 

4763350. 0 1 04 

• 

CUPP ECTI G.'lS 

1.264b 

1 . 2643 

65.4506 

aoj. cgopov> 

39074 22. 2 393 

1602396. 3662 

47 63965.4610 

V A P I A iM C E - C 0 V A R ! 

/•NCH •iATP. l.X Or The. 

STATION POSITION 



80.693*362 

-60.515006 

-4 1,4 4 4 5 33 


. -30.51300? 

94.463279 

_. 1.2., 0.7 ? 3.20 


-41.44*63$ 

12.0773 60 

33.571 770 
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« JT A~ T I ON ‘.FaR A# T= »' $ C Cfv STK *. I Nt C ' MPS7 - GSFC (GLOBAL) 



SOLUTION FOP 3TRAioL^T : !€\i.'' 1 $CAl£ *NO~ i ‘RCTATION PARAMETERS 


DX. . 

DY 

- pi 

01 CPtGA PS! 

EPSILON 

METERS 

METERS 

MET EPS 

(i 1 O'" 7 ) ' SECONDS SECONDS 

seconds' 

-10 . 6,6 

c.47 _ 

. 3 Z»;jS-Z 

9.55 0.49 0.27 

-0.40 


VAP i'aSc.E - "COVaR"! ANC? matrix’’ 


0. 6460 + 01—0.1 7 VO *01 

0.1060+01 

-0 .5130-06 

... 0.291 0-06 

0 • 1 5 « 0-06 

-0.1 76 0-07 

—0_. 17 7 Q + 01 0 . ,1,27 ' +02_ 

-0.1 6 00+ Cl 

0.1 14D-C5 

-0.10 2D- 06 

-C. 4 6 4,0-07 

-0. I 760- C6 

o.i060+oi' ; -e.;i : ss.-:+ci 

• .••o.'i.i m*cz 

— 0.6520— 06 

_0_.5 5 ID- 0.7. 

-0.126D-06 

-0.2 -5 0- 06 

_-jO . J? ,1 .S p- p 6 .0 . 1 1 4 0 - 05 

-0. c9 2p-;C6 

0.317C-12 

-0.1 120- 3 3 

-0. 5970-14 

0.2160-14 

0.29 ID- 06 r 0.1 f. 2 D- C 8_ 

.. 0..55i.D-C7 

-0 .112 0- I E 

C. 85.1 C- 13 

0.6190-14 

-0.1130-12 

0.JL.5 4Q-06.-0.46 4 0-0 7 

-o. j : eovcfc 

-C.597C— 14 

... 0.619C-1- 

0.1 04 0-12 

0. 17 3,,C-13 

0. 17?[;-(>7 -0.1 76D-C6 

-0.395D-C6 

0.2 16 C- 14 

-0.1 3 ?C- 1? 

... C.l 730-12 

0.1190-12 


CO E# I C TENT S CF CCS RE LA T I ON 


o.ioco+o: -0.1710+00 

0. 1060+00 

- 0.3 16 C +00 

0.3430+00 

0.1 640+00 

-0.1 770.-01 

■ -0.1710+00 O.iOOIT+OI 

-0. 15-.C + C0 

0.5 530+00 

-0. 176D-02 

— C . h C5 0— 01 

- 0 . i 430+00 

oi l C 6 D+ j 00 . — 0 1,1 5 40 +00 

ollOGP+iCl 

-0.3360+00 

0.551 C-UJ 

-0.1 1-0400 

-0.3VH+ )0 

-0. 3 1 6 P+Oi) 0. 56-10 +00 

- 0.35 •? C +0 0 

0.1 000+ Cl 

-C.6P2D-01 

-0.3 2 00- ci 

0.11’ 0- 01 

0.-34?n + 60?40vl7 6 r;_ C 5 

0 . 55 1 C-EC 1 

-0.6 620—01 

C. iCOl'+Cl 

o.+ssr-ci 

-0.1 IE.C + JO 

0 . 1 6 4 0 +0 p - 0 . 4 0 ? 0- C ! 

-0.11-0+00 

-0.3. VC- 01 

0.6560-0! 

O.’COD+Ol 

0.1 5-0+00 

-0.1770-01 -P.:, 4 3 0 + 00 

■•J*2 j*’ C 0 C 

D . 1 1IC-U! 

-C.l 1204 00 

c.istr+co 

C. 1 OOD+CI 





ROT AT,i ON RAP A^ETLR S CCN S 7 R 4 ! N E 13 MPS7 - SAO (GLOBAL) 


[SCt-UTriON FOR 3. TRANS l AT ION,” 1 SC ALE AND 3 FOT ATI ON PA C A ME TE PS " 


> '/:c pX" ' -OV . DZ 1. . • OL OLGA PSI E°S ILC1N 

V ■^".•••.tWETiRS . MtTfXs’t '«E-t«=‘AS‘“ ' ^,10^ St COM P S SEC ONDS SECONDS 


.. r ' * *. 

* Vi. -3 3 ' 0.52 v 25>.4i 

11.13 *:'i 

0.30 

0.53 -0. 

14 

: •. ■ 

. p'. ' : 

: ; ?? 





VARIANCE - 

COVARIANCE 

MATRIX 



o.?oio>5i 

-0.1100 + 01 0.TS7D + 00 

-0. 3*70-06 

0.1140-06 

0.1210-06 

0.1 700-07 

-o.noc+oi 

•0i?950’+Cl ••—OvSO.OD+Cl 

v ’ A ■ v • ■ : ‘ . 

0.fi?5_C-0c 

0 .1520— C 7 

-p. 2 58 0-0 7 

-0.9870- 07 

0.78 70*00 

- 0 . ?. C C D + C 1 0 . 8740+ C 1 

— 0 . 692C-.C6 

0.17COC7 

-0.6 5 CD— 07 

-0.1250-06 

— 0. 347C— 06 

}. 0.835C-C6 >09 £920-06 

0.255D -!,2 - 

0. 3800- 14 

-0.2190-14 

C.851D— 15 

0. IV, 0-06 

0 . 152D-C7 0. 13CD-C7 

-C.330C-1* 

0.3840-13 

C. 4240-14 

-0.2 67 0-14 

; 0.127 D— 06 

- 0.2 53D-C7- - C . 69 CD- G 7 

-C.219D-1A 

0 • - 2 4 C— 1 4 

0. 5 070-13 

0.5770-14 

C. 1200-S07 

. - 0 . 9 « 7 0- 07 - 0 . 1 2 5 C - 0 1 

C. 8510-15 - 

0.2870-1+ 

0. c 77D-; a 

0 . 4 6 5 P- 1 3 


» • • .s’" ! ] ‘ 






COEFFICIENTS CF CCRPFLAT! GN 









G. 1000*01 

.-Q.iji2p+oo : o;ioar+oo 

-p.259D + 00_ 

0.221-0 + 00 

0.2 030+00 

0 . 3 1 1 0- 0 1 

-0.13 2.0+00 

.0.1000+01 -0.21AD+00 

0 .5 2*0+00 

0. 2470—01 

-0. 3670-01 

-0.1*50+00 

’>•' 0.1 010+ CO 

-O.7.14D + C0 c. iocd+c: 

-C . *6 7D+ CO 

0 • w c ** i • “ 0 1 

-0.10*0+00 

“0*1 : ;7 C*-00 

rO. 2550+00 

*.0 . S2C0+q0 -0. <r'6 3C +00 

o.iooc+ot - 

0.36*0-01 

-0. 19 3D- 01 

0.7810-0? 

- 0.2? ID +00 

>0.2*70-01 0. 224D-C 1 

—0.3 8*C-C) 

O.lCOO+Cl 

OGOP-Ol 

-0. 6790-C1 

. 0.2030 *00 

-0. 3670-02 -0. 10*0+00 

-C. 1930-01 . 

0.9600—01 

O.ICOD+Ol 

0.201 n+co 

0.21;: Or 01.; 

-C . 1 *5C+CC ,-0. 177C+CC 

C.741C-C? - 

C. 6730-C1 

C.2C1D+0C 

0.130!':+ 01 
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:£ ?■ ROTATION PARAMcTERS CONSTRAINED MPS7 - GSFC (EUROPEAN) 


SOLUTION FOR 3 T R A NS L AT Yon* I SCALE AN0~3 ROTATION PAftAMETEKS 


•••• PY bz DL MEGA PS I F PS I LON 

"PETcRS METERS METERS (x 1 0" 7 ) SECONDS. SECONDS SECONDS 


-35 + 77-4 9.5 0 1 ;_4 . 5 3 -73 . 3.0_ _T 1 . 5 8 - 3 .3 5 0.03. 


• ■ ’■ • . - ■ ■ y - : - ■ ' 


' ■ :■ ' - ^ •• 

VARIANCE ■ 

- covariance 

MATRIX 



0.2340+03 

0.357D+02 

0 . 1450+03 

^a.+5.i <5 .M.Pi. 

-0.489006 

0.8570—05 

-0.1370-05 

0. 3570+02 

0.70 ID +02 

0. 9520+01 

- 0 • 5 9 0 C- 2 5 

0.6020-05 

0.2800-05 

-0.642 0-05 

0. 145D+C3 

"0. 95 20 + Cl 

0 . 2410+03 

-0.420C-04 

-0.150005 - 

■0. 1 C8D-04 

0. 2950-08 

-0.419C-04 

—0 . 5900-05 

-0.4200-04 

C. 9 3 60-1 1 

0 .1 340-12 

0.1710-12 

-C.l 590-13 

— 0. 4 P) 9 D— 06 

0.6.02D— 05 

-0. 1500-05 

0.1340-12 

0.1U'D-11 

C. 1700-12 

-0.2 350— 1 2 

; 0. 857D-05 

C. 2800-05 

-o.108r.-M: 

: 0.171 C-12 

0. 1700-12 

0.21 PD- 1 1 

— C.480D- 1 2 

-0.1870-05 

-0.6420-05 

„_0..295Dr-C5 

,C.0.„159C-i 3_ 

-0_..23f.r-12- 

C. 4 SOD-12 

0.1190-11 



COEFFICIENTS OF CORRELATION 



•• X •••. ;•••• ' fi‘ ■ 

0. 1000+01- 

0.2790+00 

0.6110+00 

— 0 .89-D + OO 

— 0 . 2 9 4 C- 0 1 

0.37SD+00 

-0.1 120+00 

0. 2790+00 

•> •• % C' ' • 

C. 100D+01 

0.73 2D- Cl 

-C. 2300+00 

0 . fc 6 2 D + 00 

0. 2260+00 

-0 .7 0? 0 + 00 

6. 61 10 + 00 

6.7320-01 

0. i.OCD+Cl' 

-0.884 O + Ou 

-0.S87D-O1 - 

o • 6 ♦ o 0 

0.1740+00 

-0. 39-C+GO 

-0.2300*00 

—0.834 n+cc 

. o.ieoc+ci 

o.4cro-oi _ 

0.3780-01 

-0.4770-0? 

-0. 2.9 4 0- Cl. 

0. 652D+00 

-0. 8 5 70- Cl 

o.405r.-c: 

0.1 OuC+01 

0.1 060+00 

-0.159C+00 

__ 1_0 -. 3 7 d 0 ♦ 0 0 : 

0.2260+00 

1 • f-/ 

-C.4o5C+0C 

0.378 C— 01 

o. ict r+co 

C. ICCD+Cl 

-0.2 970+ CO 

0. 1,1 20 ♦00 

-C. 7020+00 

o.-'r?.>r+!cd 

-0*. 4 77C-C2 

-0.1J9G+00 - 

0.2=70+00 

0. 1000+ 01 
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ROTATION PAS A METERS C CNSTR A I Nr. 0 MP37 - S AO (GLOBAL) 


S COUT TON C 0S 3 .TRANSLAT ION , 1 SCALE AND ?. R C TA f ! 0 N "p A F. 4 M E T£ RS 


vJv_ilBX __ py,- 

ME TF k S ME LE S § • ; T* £ T F R S 


-45*2 3 _.3_6 . 6A ..19 .39 


“ v VAF.LALCf - COVARIANCE MATRIX - 
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2. 5 Determination of Transformation Parameters with Constraints 

The relationship between any two geodetic reference systems would generally 
consist of seven parameters - three translations (dX, dY, dZ) between the two 
origins, three rotations (to, 0, e) of the Euler's angle type between the two sets 
of axes and the scale factor (AS), if any. 

A general transformation for the seven parameters is given below rBadekas, 

1969]: 



where to, 0 and € correspond to rotations about Z, Y and X axes respectively - the 
positive direction of rotations taken in counterclockwise mode from UVW-system to 
XYZ -system. The above equation can then be further modified as below: 



0 

1 
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0-100 
0 0-10 
10 0 -1 



0 0 -U -V 

-1 0 -V U 

0 -1 -W 0 


W 0~ 
0 -W 
-U V 




However, in the above transformation, if the geodetic reference systems are 
properly defined for Laplace condition ( parallelism of minor axis of the 
reference ellipsoid and earth's rotation axis) the three rotations arising out due to 
the improper orientation of the system are generally never more than a few 
seconds of arc while translations may amount up to 200 to 300 meters. Thus due 
to the presence of high correlation between the rotations and translations, satis- 
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factory estimates for rotations are difficult in a combined general trans- 
formation. 

An alternative method separates the determination of the rotations inde- 
pendent of the translations and the scale factor [Bursa, 1966]. The mathema- 
tical model is as follows: 

TuP- to- c cos T^tan 6^+ ^ sin T^tan 6$ = 0 

6ji^- + csinT^ + $eosT^ =0 

where T Ik and fy* are defined as the geodetic hour angle and declination of the 
(i-k)th direction of the observed point at kth station and the observer at ith station. 
The indexes (1) and (2) denote the two systems with transformation proceeding from 
system #1 to system #2. 

If we take A tt , B lk , C tt as the direction cosines of the (i-k)th direction, 

R lk as the length, then for the first system we get 



U.-P. „ Mu 

. Rik ■ . Bik ■ 



Vk...-Vi - M* 

R tt Rik 



Wy ~ W, = 6W, 

Blk Rlk 


( 4 ) 


and 



-arc tan 


Rlk 

A# 


$ik 



arc tan 





In the above relations (3,4 and 5) the elements of translati on do not enter 
the picture and a similar set of relations as per (4) and (5) above can be established 
for the second system. 
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The equation (3) then can be written as below: 
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Using the variance-covariances matrices £X andSU in respect of ith and 
kth points for the XYZ and UVW systems, the variance-covariance matrices 
S f 6 were computed for the two systems through propagation of errors as per the 
following relation [Uotila, 1967]: 

s!;> = g 

where 
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Obtaining similarly the combined variance-covariance matrix to be 
used with the equation (6) would be given by 

s!? 0 . 

0 ' ■ 

The above transformation model was used to study the relationship between 
various datums with the recent free adjustment of a Geometric Global Satellite 
Network, Solution MPS7 (fMueller, Whiting, 1972] and Section 2.4), Firstly, the 
three rotations were obtained independently of the translations with their variance- 
covariance matrices. Secondly, using the same set of common points a general 
transformation for seven parameters (including the three translations and the 
scale factor) was obtained utilizing the rotations in a constrained solution. This 
transformation was carried out in three broad groups based on the area-wise 
study i.e. , global, European and North America, with the following datums: 

(i) Goddard Space Flight Center Reference System (GSFC) r Marsh, Douglas 
and Klosko, 1971], 

(ii) Smithsonian Astrophysical Observatory's Global Reference System 
(SAO) fGaposchkin and Lambeck, 1970], 

(iii) European Datum 1950 (ED50). 

(iv) North American Datum 1927 (NAD). 

Table 1 gives the results for three rotations as obtained independently of trans- 
lations, while Table 2 gives the constrained solution for seven parameters. Table 3 
shows the results of a non-constrained general transformation for a comparative study 



The comparison shows that the constrained solutions show an overall 
improvement in all the transformations. The standard deviations in all the cases 
are smaller and the variances of unit weight show a better fit in the constrained 
solution as against the non -constrained transformation. 
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2. 6 The Impact of Computers on Surveying and Mapping 

Keynote Address Presented by Ivan I. Mueller at the Annual Meeting of the Permanent 
Committee, International Federation of Surveyors, Tel Aviv, May 29-Junc 3, 1972 

Most keynote speakers usually start with the statement that they are honored 
and privileged for the opportunity to present their views. I will not be an 
exception to this custom because I truly feel honored and privileged being 
selected by the organizing committee to deliver one of the keynote addresses 
at this meeting. Over the years, the International Federation of Surveyors has 
cortsibtfcritly sponsored a full range of valuable meetings dedicated to the exam- 
ination of important problems facing this very diversified profession. Among 
the most Innovative of the convocations called have been those associated with 
the meetings of the permanent committee. 

What then is the purpose of a keynote address ? It is generally understood to 
have a double aim. The first is to arouse unity and enthusiasm in the audience. 

But I need not concern myself with that, because I am sure that everyone here 
is equally excited at the potential of computer usage in surveying and mapping 
and at the new vistas visible on the horizon of this ancient profession. The 
other purpose of a keynote address is to present the issues inherent in the theme 
of the meeting. I shall try to present these issues, first as they are related 
to the computers, then how these machines affected traditional areas within our 
profession, what new exciting areas came into existence because the machines 
happened to be around, and finally what are those new vistas just around the 
horizon which are visible to this observer. 

The Computer 

When the computer was invented in the fifties, there was a great diversity of 
opinion on its usefulness, from skeptics who proclaimed it a toy to the more 
adventuresome prophets who predicted phenomenal growth and widespread appli- 
cation. Reflecting now on some of those early prophecies, it is obvious that 
they were vague about specific applications, real benefits , actual costs and the 
technological advances required to make the computer practical. And yet, the 


- 101 - 



usefulness has outstripped the dreams of the most adventuresome prophets. 
Undoubtedly, most people's ideas (not ours of course) about computers are 
associated with erroneous electricity or bank accounts, TV science fiction, 
moon shots or tax collection. Contrary to these beliefs, computers have a 
great deal more to offer. They work as calculators too, as repositories of 
information, as controllers, as aids to decision making in such contexts as 
banking systems, reservation systems, air and road traffic control. The use of 
computers as simulators is an application which is growing in importance: 
Examples include training astronauts, observing the effects of car crashes, 
playing war games instead of real ones, and business strategies. Computers 
have also penetrated the field of art to the dismay of some of us: Attempts have 
been made, with varying success, to use the computers as language translators, 
as writers of poetry and prose, as producers of visual art, to create ballet 
routines, and both write and synthesize music. There is plenty of scope here 
for those of us who enjoy a debate guaranteed to have no conclusive outcome. 

On the serious side, because of its varied applications, the computer demands 
from society, including the surveyors, decisions as important as any it has made, 
certainly as important as those forced on our predecessors by the industrial 
revolution. It is sad that the level of discussion, even in some "professional" 
circles, has so far been so puerile, to understanding of the issues so limited 
and so inadequate. 

With this in mind, allow me, in a few minutes, to review the progress over the 
past two decades to see how the use of computers has developed and then to 
examine current trends. 

The first decade of computer development, in the 1950's, saw the use of 
machinery largely as an aid to scientific research; many research projects in 
physics, chemistry and engineering demand elaborate calculations - the design 
of an aircraft wing or engine, for instance, or the design of a nuclear reactor. 

As a matter of fact, there is one project - atomic bomb development - which has 
always demanded more and more calculations in order to progress with as little 
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testing as possible. It is easier and also rather more socially acceptable to 
simulate an explosion on a computer, however large and expensive, than to 
explode a live bomb. This one use played an important part in the development 
of very large and very fast computers during that first decade. It was not until 
after 1960 that such machines found their way into other than atomic research 
laboratories. The second decade of computer development, in the 1960's, 
saw the development of the computer as an electronic office, a data handler and 
processor. The computers initially used in this era were designed not as 
calculating engines for scientific use but to make the processing of card files 
cheaper and easier. The jobs being done were those which are carried out 
within the administrative and accounting departments of a business. Such jobs 
placed more emphasis on the storage capacity available in the machine than on 
its calculating speed - in contrast to the research applications in the first decade. 

As the users became more confident in and more used to computers, new 
applications appeared using both the calculating capacity of the machinery and 
its data handling capabilities. 

In looking back, it becomes relatively easy to separate the demarcation points 
between post generations of computers. Historically, these have occurred following 
advances in hardware technology: vacuum tubes for the first generation around 
1950/51, transistors for the second (between 1958 - 60), and integrated transistor 
circuits for the third between 1963 and 1965. Lately, however, the introduction 
of many other new features - in peripherals, communications, remote terminals, 
operating systems, and the like - have made the distinction between the generations 
increasingly fuzzy. We have now passed the eve of the fourth generation computers 
which is best characterized by the ability to provide information which is constantly 
on the tap. In other words, while the roles of the first three generations were 
computations, data and information processing, the current generation also pro- 
vides on-line information. The rapid evolution through the fourth generation - 
spurred on primarily by the immense proliferation of minicomputers - is under- 
way and one can now begin to imagine the hardware and software components 
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which will characterize the fifth generation projected to be bom between 1975 
and 1978. 

I will not elaborate on the technical aspects of these future babies of the 
computer industry. Let me just say that these new machines are being viewed 
as man’s "intelligent” assistants. Many of them will be portable, hand carried 
or in the car and in the home, that can be plugged into telephone and electric 
outlets or even carry their own power supply. This will tie the computer 
completely to the telecommunications systems, allowing the computer to 
’remote' its power to where it is needed. Indeed, the telephone will become 
probably the most widely used terminal of the 1970's - incorporating voice 
output and touch tone input. Such an availability of computer power can have 
nothing less than an immense impact on society, greater perhaps even than 
the impact television has had. 

New major innovations are likely to occur also in the software area. For 
instance,the cost of programming, which has been held almost constant (per 
line of code) throughout the past three generations, should be reduced by more 
than a factor of ten in fourth generation systems. This should come as a direct 
result of interactive programming using time shared facilities. A further factor 
of ten reduction in costs can be expected with the fifth generation. With the 
remote terminal and the packaged programs (to which I will return a little later) 
will come a truly conversational use of the computers. Many such systems are 
now being designed and use languages suitable even for the non-professional. By 
the end of the fifth generation - by the early 1980's - literally anyone will be able 
to use a computer and many programs should be available for helping us perform 
our daily tasks. Computers and terminals could then become as common as 
telephone and television today. 

In passing through the second and third generations of computers there was 
approximately a four fold increase in the number of computers in use per generation. 
Throughout the 1960's there was a ten fold increase. Assuming that these trends 
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continue, then by 1975 - at the onset of the fifth generation - there will be more 
than 200, 000 computers in use around the world. By 1980 there could be over 
500, 000. But if we count the remote terminals, then these numbers grow by a 
further factor of ten. Moreover, if we include all the telephones used for 
remote access to computers, then practically everyone with a telephone will 
have access to a computer by 1980. 

What are the uses of all these computers? In addition to applications in our 
own profession there are of course countless applications. Let me select for 
illustration probably the most sophisticated one, the applications in management 
science: 

The major object of modern computer applications in this field is the setting 
up of a computerized data base to enable better analysis to be made of alternative 
uses of resources. At present, many important decisions are taken on inadequate 
data or on information which is out of date. In a stable and well-established 
business this may be of little consequence, but for firms in rapidly changing 
markets or involved in rapid growth or technological change, timeliness of data 
can be vital. Rapid and convenient access to the data base is therefore required, 
and it is necessary that the whole system be designed so it can react to the users 
urgent demands. Modern computer techniques enable the user to converse with 
the computer over a terminal. The user can ask questions of the computer, 
which can then, by questioning the user, ellicit further information to retrieve 
the answers required from its memory. In this way, the data base can be 
searched, and the result of a requested analysis can be made instantaneously 
available. 

The nature and complexity of the analysis required may differ considerably, 
go that it would be inefficient to have the most powerful processor tied up wholly 
with one user. The equipment needed to implement such an enquiry system is thus, 
not one computer, but a collection of units, some of which are devoted mainly to 
manipulating data, some to the calculations needed for analysis of the data, some 
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to up-dating the files as new data arrives, and some to conversing with the 
users. As the users and the data sources may be physically distributed over 
a wide geographic area, the whole complex must be connected by communications 
channels, and thus, becomes a computer-network. At present such networks 
are being built for several applications. Several already exist - for instance, 
to carry out airline and hotel reservations on a world-wide basis. Others are 
being installed to link hospitals into the data base containing information on 
patients, availability of beds, etc. There is no intrinsic reason why, in due 
course, single overall systems should not serve the needs of all the users in 
any technical or geographical group desired. Several computer bureau operators 
with machines in different countries are planning to link their machinery so 
they can work on whichever machine is most readily available or most economic 
at the time. Such arrangements could well form the basis for an international 
computer-network.. 

The establishment of such a network naturally will contain some inherent 
dangers for the individual, primarily related to his status within the community, 
who can be affected without his knowledge. In order to bring about beneficial 
applications, the computer must have data - not only about money and materials 
and the rest of the physical environment in which we live, but also about people 
and their attitudes and circumstances. Until recently the clerical effort needed 
to cross-reference all these files has fortunately been prohibitive. But once 
these data find their way into a computer system, cross-connections could be . 
made in a matter of second. Thus, on applying for an insurance benefit you 
might find the amount of your last unpaid parking fine deducted automatically, 
or perhaps find yourself arrested to answer a charge of speeding. Would we 
be happy under an efficient tyranny - one in which every movement and action 
of the citizen was recorded, analyzed, cross-checked instantaneously and no 
incident, no matter how trivial, ever forgotten? Yet, such is the mechanism we 
now have the capacity to create. It is not a far stretch of the imagination from 
here to see that Orwell’s 1984 predictions on surveillance could also be fulfilled 
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and on schedule. 

It is not, of course, the computer itself which creates social problems, 
but the human beings into whose hands it is placed. The computer is a tool 
and it can be used or abused by man at his discretion. Compared with such 
tools as nuclear energy, the computer does seem to possess more potential for 
good than harm. 

Whether this picture appeals to you or frightens you, I have no way of knowing. 
A recently published book entitled Future Shock, concerned itself with the plight 
of modern man attempting to cope with "an environment so ephemeral, unfamiliar 
and complex as to threaten millions with adaptive breakdown. " The book is an 
indication of the apprehension with which some people view the future and it is 
worthwhile for those of us who are contributing agents of technological evolution 
to do some hard thinking about where we are going - to alleviate the fears of 
some and help all prepare for the coming advance in technologies. 

Let us now take a look at how the availability of the generations of computers 
affected surveying and mapping. Obviously, this review will have to be a selective 
and a subjective one. I will be able to mention only the most spectacular examples 
and only those which are likely to be in the interest of this convocation, and of 
course, only those which are in my area of competence. 

The Shape of the Earth and its Gravity Field 

I should make it clear at the outset, that I am not concerned with local ir- 
regularities in the earth's surface, the mountains and the valleys. I shall be 
discussing the mean sea level surface of the earth, carried through under the 
land, the surface usually called the geoid. This geoid, being a surface on which 
the potential of the earth's gravity field is constant, will, at the same time serve 
as a pictorial representation of the variations in the gravity field of the earth as 
well. 

In the United States, a historical review on the subject in "which shape the 
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earth is in, " probably would start from the time when the Declaration of 
Independence from England was signed in 1776. In this country, however, one 
is obliged to start with the prehistoric man, who, if he thought about the subject 
at all, presumably concluded that, apart from local oddities like rocks or 
mountains, the earth was flat. This is also the view held today by the Flat-Earth 
Society, also in England. 

The idea of. a nearly-spherical earth was surprisingly late in becoming 
established, or so it seems to us, with the advantage of hindsight. Neither the 
Babylonians nor the Egyptians favored dlls idea, and the credit goes to Pythagoras 
and his school in the sixth century B.C. I should add that the Idea was derived not 
from observations but from their conviction that the sphere was "the perfect" 
shape. Three-hundred years later Eratosthenes did more than adopt the idea, 
he actually measured the earth's circumference, using the propagation velocity 
of a camel caravan as his scale. 

It was not until the seventeenth century that the shape of the earth was improved 
upon. The first indication that the earth may be flattened at the poles was obtained 
in 1672 by Jean Richer’s French expedition to South America, where he found 
that his pendulum clock, accurate in Paris, was loosing time at Cayenne. First 
numerical estimates on the flattening came from Newton in his "Principia" 
published in 1687, but practical measurements to establish the value of the 
flattening were made by the Cassinis, who measured arc length in France, 
and who came to the conclusions that the earth was flattened indeed, but not at 
the poles, but rather at the equator, thus, it looked like an egg or a lemon. This 
was in 1720, and . fierce controversy followed: Was the earth flattened or elon- 
gated at the poles? Who was right, Newton or the Cassinis? The French 
Academy sent the two famous expeditions of Maupertuis to Lapland, and La 
Condamine's to Peru. After ten years of labor and an equal number of years 
spent in quarrels, the conclusions tended to confirm Newton's idea, and 
Voltaire congratulated the expeditions saying, "You have successfully flattened 
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the poles and the Cassinis. " Not much happened after the regarding the shape of 
the earth, until the middle of the present century, when first analysing gravity 
measurements on a global basis, and after 1957, analysing the orbits of 
artificial satellites, a complete new picture of the earth's shape emerged. 

These analyses, of course, were made possible only because by that time, the 
computers came into existence. 

I shall not describe how from the perturbations to satellite orbits, caused 
by the various possible oddities in the earth's shape, these oddities can be deter- 
mined. It should suffice to say that a new value for the flattening has emerged, 
indicating that the earth's equatorial diameter exceeds the polar diameter by 
42. 77 km, which is a full 170 meters different from the previously adopted 
value. This difference may not seem much for most of us, but it is important 
for the geophysicist, who may conclude that the earth's interior has great 
strength, and the assumption that it can be treated as if it were a fluid, an 
assumption which in the past, was. widely made, is illegitimate. 

The more accurate value for the flattening is, however, only a very small 
part of the information obtained from satellites. Without going into technicalities, 
let me simply illustrate the improvement by the fact, that in the pre-satellite era, 
the shape of the earth and its gravity field was described by four basic parameters, 
while today, the number of known parameters exceeds two-hundred and fifty. This 
new information pictorially represented as the aforementioned geoid above the 
ellipsoid shows that the most prevailing features are the healthy depression 
around the South Pole, a bulge south of the equator, and also around the North Pole, 
indicating, in the language of the press, that the earth is "pear-shaped. " This 
discovery came as a relatively great surprise to most of us, but it should have 
been no surprise to Christopher Columbus, who gave it as his opinion "that 
it has the shape of a pear that is very round, except where the stem is, which is 
higher. . . " Other important features are the depression south of India, (113 m), 
the elevation near New Guinea, (81m), and the elevation centered in England 
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and the south Atlantic. 

To sum up, satellites and the computers have brought us from the earth of 
1957, which was merely a sphere flattened at the poles, and flattened by the 
wrong amount, to a complicated figure which when seen in the round looks 
perhaps like a potato with dips and humps all over it. 

By-products of this satellite-orbit analysis are the coordinates of the tracking 
stations with respect to the center of the earth. In the pre-satellite era, such 
information, which is vital in relating the numerous geodetic systems of the 
world, practically did not exist. Today, geocentric coordinates are known 
for about 150 stations fairly evenly distributed around the globe. 

Satellites also help in mapping, as geometric triangulation points in the sky 
in connection with the method called: 

Satellite Triangulation or Trilateration 

This method found wide range applications in connecting another 150- 200 
tracking stations in the relative sense both on a continental and on a global 
basis. Better known projects in this category are the programs under the coor- 
dination of the Eastern and Western European Subcommissions for Satellite 
Triangulation of the International Association of Geodesy; the U.S. National 
Geodetic Satellite Program now in its final stages, including observations by 
the Smithsonian Astrophysical Observatory, The National Geodetic Survey, 
(formerly Coast and Geodetic Survey), NASA and various other agencies; the 
French coordinated ISAGEX Program; other French works in southern Europe 
and northern Africa; and some other local national network developments in 
North and South America. 

I will not attempt to offer you a glimpse at the software used in the calculations 
related to satellite geodesy, mainly because some of them are rather lengthy. 

The fact that some of these programs took 100 man years to develop is an indi- 
cation not only of the complexity of the problem, but also of the need for better 



programming methods. Clearly, when one needs to work with several 
ten-thousand observations in order to determine several hundred unknown 
quantities, like station coordinates, gravitational parameters, and at the same 
time, attempts to recover at least some of the systematic errors burdening 
the observations, the computer software and hardware will have to be impressive 
indeed. 

This leads us to an application where the impact of computers is and will 
probably be the greatest both in its economical aspects and also in the number 
of people affected. This application is generally known in surveying circles as: 

Adjustment Computations 

Adjustment in the surveying and mapping terminology is the method used to 
derive unique and "best” values for parameters from redundant measurements 
of those parameters, or parameters related to them by a known mathematical 
relationship. It is a device which should be used by everyone in the profession 
involved in the evaluation of survey data from leveling to satellite laser ranging 
or, from cadaster surveys to lunar mapping. The fundamentals of this science 
were laid down by Karl Friedrich Gauss in the eighteenth century at the age of 18. 
Every geodesist and photogrammetrist of note since then, has contributed to 
the literature by refining (or confusing) some aspects of the topic. 

Without going again into the technical details to the extent possible, let me 
remind you that in the pre-computer era, up to the early fifties, one did not 
enter lightly into an adjustment computation; one looked very closely at the 
model; one checked and double-checked the input data, and in very special cir- 
cumstances, one might undertake the extra computations necessary to check the 
possible correlations between the unknown parameters, or to compute the error 
ellipses for certain selected points of special interest. In other words, it was 
not practically feasible to put adjustment computations on a sound statistical basis. 
The number of unknown parameters was also limited, since the computations 
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for one medium-sized network (50 - 100 unknowns) were likely to require 
several man-months of time, thus, it was a very expensive undertaking indeed. 
The use of statistical methods for planning a network to make sure that it is 
the most economical and most favorable from the point of view of the propagation 
of errors was almost out of the question because of the costs Involved. For 
this reason, in a given country, very few organizations were doing adjustment 
computations. 

Today, thanks to the computers, this situation is part of history. Very large 
numbers of organizations are doing adjustment computations using computer 
programs, either developed by themselves or procured from other organiza- 
tions. These programs are (or should be) based on correct statistical theory 
and techniques, and running them, even with a very large number of unknowns, 
costs very little. 

Advances in this regard were most spectacular in that part of the mapping 
industry which deals in photogrammetry, where the wide applications of aerial 
triangulation or analytical photogrammetry using block adjustment techniques 
with a great number of unknowns is part of the daily routine. Another spectacular 
area where adjustment computations are routinely used to full capacity is 
satellite geodesy, where the number of unknown parameters, mostly highly 
correlated, and to be adjusted for in one huge simultaneous adjustment, may 
reach several thousand. 

It is interesting to note that a significant number of rather sophisticated 
"package programs" written for different purposes, like aerial triangulation, 
horizontal control, satellite triangulation or orbit determination, have been widely, 
distributed and used by a great number of organizations other than those who 
designed the programs. It is a small step from here to arrive to the point, where 
the average surveyor can pick up his phone and dial the computer or go to his 
remote terminal, specify his object, read the input data in the specified manner, 
and receive his results with all the statistical trimmings faster and cheaper 



than he ever dreamed of. He has a powerful design tool at his command; he 
can now make full use of law of error propagation and optimize any system he 
is designing; he can build in constraints; he can test options and find the 
option that meets his specifications with the least effort and cost. At the con- 
clusion of the project, he can do an evaluation and test the assumptions that it 
was necessary to make about his instruments. If data from a variety of sensors 
have been combined in an adjustment, he can test the distribution of residuals 
for normality; he can test his mathematical model, his weighing procedure. 

In theory, this always has been possible, but until modem computer facilities 
became available, it was out of the question as a regular tool. 

Equipment Oriented Areas 

There are also equipment oriented areas where the availability of the gen- 
erations of computers (directly part of, or tied to, the sensor-system) affected 
surveying and mapping. To mention a few, let me start with the AN/USQ-28 
Mapping and Surveying System, which comprises the most advanced group of 
equipment integrated to collect accurate raw data for mapping purposes. It 
was specifically designed to acquire photography suitable for 1:50, 000 scale 
topographic mapping in areas where ground control is insufficient. The system 
is built into a Boeing 707 aircraft and consists of precision mapping cameras, 
an inertial navigation system, electronic distance measuring equipment, a terrain 
profile recorder, and other auxiliary equipment. All data, with the exception of 
the photography, are recorded on magnetic tapes for direct input into digital equip- 
ment to speed the data reduction process. It is a pity that as of this moment, 
the system is not operating because there seems to be lack of money to pay for 
the operation of the aircraft (for gasoline!). 

Another example is the progress that has been made in automated computation 
equipment. These computer-driven machines use image sensing and correlation 
techniques to produce horizontally correct images while simultaneously detecting 
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and recording height information. This equipment is supposed to reduce map 
compilation time by 75%. 

Another development of significant interest is the automatic or semi- 
automatic coordinate readers. This equipment is designed to measure, for 
example, the coordinates of star images on photographic plates obtained for 
astronomic or satellite geodetic applications. The instruments have a pre- 
programming feature which moves a detecting head tQ the approximate location 
of each required star image. The detection head then centers itself precisely 
over a star image, at which point the coordinates are measured and recorded 
on punch cards for input in the computer program. 

Another and rather esoteric computer application in this equipment oriented 
category is the Apollo mapping system for accurate lunar mapping. The main 
purpose of the system is similar to that of the USQ-28 mapping system mentioned 
earlier, i.e. , to provide maps in areas where ground control is insufficient. 

The lunar orbiter and Apollo programs through Apollo 14 have produced phenom- 
enal photography to support landing site selection and surface operations. How- 
ever, the new metric camera system which was flown first on Apollo 15, then 
on Apollo 16 and which will be also on board the last manned flight to the moon, 
offers an order of magnitude improvement towards lunar mapping, the deter- 
mination of the lunar gravity field, and of the motion of the moon in space. It 
is again a pity that the system is included only in the last three missions, and 
was left out from the previous seven missions. Of course, the astronauts on 
Apollo 7-10 were rather busy preparing the landing of Neal Armstrong onApollo 
11, but only NASA knows why the system was not flown on Apollo 12-14. The 
area coverage would have been certainly better. 

This system consists of three cameras, a laser altimeter and timing equipment. 
The first camera is a 3-inch metric mapping camera which photographs the lunar 
surface while the second stellar camera built into the same housing takes simul- 
taneous pictures of the star field just above the lunar horizon to aid the deter- 
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mination of the orientation of the mapping camera. The laser altimeter is 
synchronized to fire simultaneously and provides the distance from the camera 
to the lunar ground for each photograph. All this information together with 
the earth-based tracking data should give sufficient information on the position 
and orientation of the mapping camera (to about 2.5 m relative). The third 
24 inch panoramic camera provides very high resolution photographs (2 m at 
the nominal 110 km altitude). 

As I mentioned, the main application of the system is to establish geodetic 
control on the moon and provide maps for the areas covered. In addition to 
these, information is expected on the rotation of the moon about its axis, com- 
monly known as the phenomena of libration. The data will also be analyzed in 
conjunction with the laser distances measured between earth-based observatories 
and the reflectors placed on the moon surface by Apollos 11, 14, and 15, and 
Luna 17. This combination of data should be most helpful to improve on the lunar 
ephemeris, i. e. , on the knowledge of the relative motion of the moon around the 
earth, which lately seems to be part of geodesy also. 

The Future 

From here, there is only a short step into the future. What will the next 
decade bring? I already described what is expected from the computers and how 
they will change the job of the surveyor in the adjustment area. Let us see 
briefly that in addition to the routine mapping and surveying activities, what 
miracles the surveyor is to perform during the next decade or so. First of all, 
he is going to get some new customers: the geophysicists and the oceanographers. 

He will need new tools, because their demand for a full magnitude and better 
positions (from 10m to I'm to 10cm) than what is available today exceeds present 
capabilities. Most of these instruments are already in the development stage 
and undoubtedly will be ready for applications in the not too distant future. Let 
us take a quick look at these machines: 
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First of all, existing laser distance measuring devices will be improved to 
the point where the only factor limiting the accuracy of the observations will 
be the uncertainties in tropospheric propagation, which is expected to be reduced 
to about 6 cm (from the present 15 - 30 cm). 

On the radio frequency systems with prospects of 1 m or better accuracy, 

Very Long Base-Line Interferometry (VLBI) seems to offer the greatest 
versatility. This technique depends upon local frequency standards of high 
quality - preferably hydrogen masers - at two or more radio antennae separated 
by distances on the earth as great as allowed by the common visibility of a radio 
source, like a quasar or a water vapor source. The frequency standards provide 
time references for magnetic tape recordings of signals from these galactic 
energy sources. The tapes are later correlated at a central computing facility, 
and the time difference for arrival of the same wave is determined. From this, 
it will be possible to calculate the distance between the two antennae to an 
accuracy of about 15 cm and the direction between them to about 0. 001 arc second, 
provided that the position of the energy source Is known. 

Satellite to satellite (range rate) tracking also will offer substantial advantages 
over current techniques limited by our dirty window towards space, the atmos- 
phere. Very high satellites will track a low satellite continuously through the 
vacuum of the universe with very high precision. 

Such continuous tracking technique, coupled with the so-called "drag-free" 
satellite, will further improve our knowledge of the gravity field of the earth and 
the geoid. The essential element of such a system is an unsupported mass 
contained in a spherical shell. A control system in the satellite senses motions 
of the shell relative to the proof mass and actuates small thrusters that force 
the shell to follow the proof mass without touching it. Hence, the proof mass 
follows an orbit influenced only by gravitational force. 

Improvement in the knowledge of the gravity field, the shape of the geoid is 
also expected through the satellite to ocean radar altimeters, measuring contin- 
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uously the distance between sea level and the satellite. The first of these devices 
will be flown probably in 1974 on an experimental basis. 

From these new instruments, a wealth of information will be made available 
to the earth scientists, who, in turn, will be able to produce unpredictable 
but certainly substantial advances about the rotational motions of the earth, 
tide interactions, temporal variations in the gravity field, continental drift and 
other large scale deformations of the earth crust and mantle. The interactions 
of these motions and deformations appear to be responsible for a wide variety of 
effects, including large earthquakes, mountain building, generation of tsunamis 
(tidal waves) , and confinement of nearly all active volcanoes to only a few narrow 
belts. The satellite born radar altimeter eventually will provide valuable ocean- 
ographic information on tides, storm surges, general ocean circulation, and 
other dynamical processes affecting sea level. 

Most of these problems are global in nature, thus, require observations 
globally distributed. The interaction between the several dynamic subsystems of 
the earth demands coordination of the observations. Hence, for maximum effec- 
tiveness, technological integration and international cooperation are essential to 
a progressive investigation of these topics. 

Is the International Federation of Surveyors willing and ready to participate in 
this cooperation? It it ready and willing to take this challenge and serve the new 
customers ? 

What else is coming? - Automated data banks with national and international links 

- Automated data reduction systems 

- Remote sensing satellites for environmental monitoring, 
ocean sensing and for land use and resources management, producing 15, 000, 000 
bits of information per second - equivalent to an Encyclopedia Britannica every 
couple of minutes. We certainly will be able to verify the conjuncture that as civi- 
lized man evolved from his primitive ancestry, he developed an appetite for large 
masses of data, recording observations about his individual or collective activities 
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with ever greater precision and detail. 


Concluding Remarks 

On the surface, it would seem that surveying presents no serious issues as 
a technology. It is a useful tool in the service of mankind and extends the capa- 
bilities of science. Unlike some technologies , surveying does not pollute. On 
the contrary, it may help to preserve the quality of the environment. It is not 
likely to be wasteful economically. Instead, it could stimulate and guide resource 
development as well as scientific research in the earth sciences. What is more, 
it has some popular attributes. It requires a private and public sector team 
effort, and is multi-disciplinary as well as multi-institutional and multi-national. 

But, if we were to conclude from such reasoning that no major issues are 
involved, we would be badly mistaken. The issues are not technological, but 
sociological. In my view, they effect the unity of the profession of surveying 
and mapping. 

Let me quote a recent editorial from the transactions of the AGU on the 
"Surveyor Geodesist": 

"For over two thousand years, the land surveyor and the astronomer 
often joined by the mathematician, collaborated in the development 
of geodetic science. This symbiotic relationship, which reached its 
zenith in the last three hundred years, resulted in inferences of geodetic 
significance from observational data and also led to the establishment 
of the science on a rigorous mathematical foundation. The surveyor, 
to some degree and to a limited accuracy, participated in this develop- 
ment in the small; but, today he is severely hampered by the restrictive 
limits to his data base, by the limited scope of his observing instruments 
and computing methods, and, in no small way, by the deemphasis in 
surveying education at the university level. In addition, photogram- 
metric methods and, in more recent times, developments in space 
technology have made enormous inroads into his areas of competence. 

In fact, the phenomenal geodetic fallout from the space program has so 
obscured the place of the surveyor in the geodetic scheme of things that 
there is a tendency to downgrade his continued vital contribution to the 
science. Hence, more and more the average surveyor finds himself 
outside the geodetic mainstream, relegated to a supporting role as a 
provider of cadastral and lower order engineering data. 



The new team combines the expertise of the mathematician, the 
physicist, and the space scientist. From space-oriented observations, 
this group of scientists has obtained data in regions inaccessible to 
the surveyor and has obtained results that the geodesist using classical 
techniques could never hope to achieve. As the space scientist refines 
his measurements and increases his sampling rate, thereby providing 
more precise data at ever decreasing wave lengths, the geodesist finds 
that among many applications he can support the oceanographer in 
resolving ocean surface problems; the tectono-physicist and the 
seismologist in measuring continental drift and crustal movement; 
and the astronomer in determining polar motion and variations in earth 
rotation. 

Will this expanded geodetic role further divorce the surveyor from 
the geodetic community ? Not necessarily; a great deal depends upon 
the willingness of the profession to broaden its horizons. The new 
users of geodetic information require baseline information at accuracies 
comparable to and sometimes exceeding those the surveyor is accustomed 
to providing on a day-to-day basis. The surveyor needs to seek out his 
new customers and needs to become aware of his problems; he needs to 
upgrade his field operations, using the most precise instrumentation 
and adjustment techniques; and he most certainly must insist upon 
improving and expanding the university curriculum in surveying." 



3. PERSONNEL 


Ivan I. Mueller, Project Supervisor, part time 
Muneendra Kumar, Graduate Research Associate, part time 
James P. Reilly, Graduate Research Associate, part time 
Narendra K. Saxena, Research Associate, full time 
Tomas Soler, Graduate Research Associate, part time 
Emmanuel Tsimis, Graduate Research Associate, part time 
Marvin C. Whiting, Gradute Research Associate, part time 
Susan Breslow, Research Aide, part time 
Barbara Beer, Research Aide, part time 
Evelyn E. Rist, Technical Assistant, full time 


4. TRAVEL 

Ivan I. Mueller 

Bronx, New York, February 22, 1972 

Attend GEOP Research Conference Steering Committee Meeting 
Ivan I. Mueller 

Graz, Austria, May 29-June 2, 1972 

Present a paper at IAG International Symposium on Satellite and 
Terrestrial Triangulation 



5. REPORTS PUBLISHED TO DATE 

OSU Department of Geodetic Science Reports published under Grant 
No. NSR 36-008-003: 

70 The Determination and Distribution of Precise Time 
by Hans D. Preuss 

April, 1966 

71 Proposed Optical Network for the National Geodetic Satellite Program 
by Ivan I. Mueller 

May, 1966 

82 Preprocessing Optical Satellite Observations 
by Frank D. Hotter 
April, 1967 

86 Least Squares Adjustment of Satellite Observations for Simultaneous 
Directions or Ranges, Part 1 of 3: Formulation of Equations 

by Edward J. Krakiwsky and Allen J. Pope 
September, 1967 

87 Least Squares Adjustment of Satellite Observations for Simultaneous 
Directions or Ranges, Part 2 of 3: Computer Programs 

by Edward J. Krakiwsky, George Blaha, Jack M. Ferrier 
August, 1968 

88 Least Squares Adjustment of Satellite Observations for Simultaneous 
Directions or Ranges, Part 3 of 3: Subroutines 

by Edward J. Krakiwsky, Jack Ferrier, James P. Reilly 
December, 1967 

93 Data Analysis in Connection with the National Geodetic Satellite Program 
by Ivan I. Mueller 
November, 1967 

OSU Department of Geodetic Science Reports published under Grant 
No. NGR 36-008-093: 

100 Preprocessing Electronic Satellite Observations 
by Joseph Gross 
March, 1968 

106 Comparison of Astrometric and Photogrammetric Plate Reduction Techniques 
for a Wild BC-4 Camera 
by Daniel H. Hornbarger 
March, 1968 


- 123 - 



110 Investigations into the Utilization of Passive Satellite Observational Data 
by James P. Veach 
June, 1968 

114 Sequential Least Squares Adjustment of Satellite Triangulation and 
Trilateration in Combination with Terrestrial Data 
by Edward J. Krakiwsky 
October, 1968 

118 The Use of Short Arc Orbital Constraints in the Adjustment of Geodetic 
Satellite Data 
by Charles R. Schwarz 
December, 1968 

125 The North American Datum in View of GEOS I Observations 
by Ivan I. Mueller, James P. Reilly, Charles R. Schwarz 
June, 1969 

139 Analysis of Latitude Observations for Crustal Movements 
byM.G. Arur 

June, 1970 

140 SECOR Observations in the Pacific 

by Ivan I. Mueller, James P. Reilly, Charles R. Schwarz, Georges Blaha 
August, 1970 

147 Gravity Field Refinement by Satellite to Satellite Doppler Tracking 
by Charles R. Schwarz 

December, 1970 

148 Inner Adjustment Constraints with Emphasis on Range Observations 
by Georges Blaha 

January, 1971 

150 -Investigations of Critical Configurations for Fundamental Range Networks 
by Georges Blaha 
March, 1971 

177 Improvement of a Geodetic Triangulation through Control -Points 
Established by Means of Satellite or Precision Traversing 
by Narendra K. Saxena 
In press 


- 124 - 



The following papers were presented at various professional meetings: 

’’Report on OSU participation in the NGSP" 

47th Annual meeting of the AGU, Washington, D. C. , April 1966 

"Preprocessing Optical Satellite Observational Data" 

3rd Meeting of the Western European Satellite Subcommission, LAG, Venice, 
Italy, May 1967. 

"Global Satellite Triangulation and Trilateration" 

XIVth General Assembly of the IUGG, Lucerne, Switzerland, September 1967, 
(Bulletin Geodesique , March 1968). 

"Investigations in Connection with the Geometric Analysis of Geodetic Satellite 
Data" 

GEOS Program Review Meeting, Washington, D. C. , Dec. 1967. 

"Comparison of Photogrammetric and Astrometric Data Reduction Results for 
the Wild BC-4 Camera" 

Conference on Photographic Astrometric Technique, Tampa, Fla. , March 1968. 

"Geodetic Utilization of Satellite Photography" 

7th National Fall Meeting, AGU, San Francisco, Cal. , Dec. 1968. 

"Analyzing Passive-Satellite Photography for Geodetic Applications" 

4th Meeting of the Western European Satellite Subcommission, LAG, Paris, 

Feb. 1969. 

"Sequential Least Squares Adjustment of Satellite Trilateration" 

50th Annual Meeting of the AGU, Washington, D. C. , April 1969. 

"The North American Datum in View of GEOS-I Observations" 

8th National Fall Meeting of the AGU, San Francisco, Cal. , Dec. 1969 and 
GEOS-2 Review Meeting, Greenbelt, Md. , June 1970 (Bulletin Geodesique, 
June 1970). 

"Experiments with SECOR Observations on GEOS-I" 

GEOS-2 Review Meeting, Greenbelt, Md. , June 1970. 

"Experiments with Wild BC-4 Photographic Plates" 

GEOS-2 Review Meeting, Greenbelt, Md. , June 1970. 

"Experiments with the Use of Orbital Constraints in the Case of Satellite Trails 
on Wild BC-4 Photographic Plates” 

GEOS-2 Review Meeting, Greenbelt, Md. , June 1970. 


- 125 - 


"GEOS-I SE COR Observations in the Pacific (Solution SP-7)" 

National Fall Meeting of the American Geophysical Union, San Francisco, 
California, December 7-10, 1970. 

'Investigations of Critical Configurations for Fundamental Range Networks" 
Symposium on the Use of Artificial Satellites for Geodesy, Washington, D. C., 

April 15-17, 1971. 

"Gravity Field Refinement by Satellite to Satellite Doppler Tracking" 

Symposium on the Use of Artificial Satellites for Geodesy, Washington, D. C. , 

April 15-17, 1971. 

"GEOS-I SECOR Observations in the Pacific (Solution SP-7)" 

Symposium on the Use of Artificial Satellites for Geodesy, Washington, D. C. , 

April 15-17, 1971. 

"Separating the Secular Motion of the Pole from Continental Drift - Where and 
What to Observe?" 

IAU Symposium No. 48, "Rotation of the Earth, " Morioka, Japan, May 9-15, 1971. 

"Geodetic Satellite Observations in North America (Solution NA-8)" 

Annual Fall Meeting of the American Geophysical Union, San Francisco, 
California, December 6-9, 1971. 

"Scaling the SAO-69 Geometric Solution with C-Band Radar Data (Solution SC 11)" 
Annual Fall Meeting of the American Geophysical Union, San Francisco, 
California, December 6-9, 1971. 

"The Impact of Computers on Surveying and Mapping" 

Annual Meeting of the Permanent Committee, International Federation of Surveyors, 
Tel Aviv, Israel, May 1972. 

"Investigations on a Possible Improvement of Terrestrial Triangulation by Means 
of Super-Control Points" 

IAG International Symposium - Satellite and Terrestrial Triangulation, 

Graz, Austria, June, 1972. 

"Free Adjustment of a Geometric Global Satellite Network (Solution MPS7)" 

IAG International Symposium - Satellite and Terrestrial Triangulation, 

Graz, Austria, June, 1972. 


- 126 - 


